TASI Lectures on Inflation 



William H. Kinney 
Department of Physics, University at Buffalo, SUNY, 
Buffalo, NY 14260-1500, USJ^ 

This series of lectures gives a pedagogical review of the subject of cosmological inflation. I discuss 
Friedmann-Robertson- Walker cosmology and the horizon and flatness problems of the standard hot 
Big Bang, and introduce inflation as a solution to those problems, focusing on the simple scenario 
of inflation from a single scalar fleld. I discuss quantum modes in inflation and the generation of 
primordial tensor and scalar fluctuations. Finally, I provide comparison of inflationary models to 
the WMAP satellite measurement of the Cosmic Microwave Background, and briefly discuss future 
directions for inflationary physics. The majority of the lectures should be accessible to advanced 
undergraduates or beginning graduate students with only a background in Special Relativity, al- 
though familiarity with General Relativity and quantum field theory will be helpful for the more 
technical sections. 



I. INTRODUCTION 



Cosmology today is a vibrant scientific enterprise. New precision measurements are revealing a universe with 
surprising and unexpected properties, in particular the Dark Matter and Dark Energy which are now believed to 
be the dominant components of the cosmos. Galaxy surveys such as the Sloan Digital Sky Survey are making the 
first large-scale maps of the universe, and satellites such as WMAP are making exquisitely precise measurements of 
the Cosmic Microwave Background (CMB), the haze of relic photons left over from the Big Bang. In turn, these 
measurements are giving us clues which are helping to unravel one of the oldest and most profound questions people 
have ever asked: Where did the universe come from? In these lectures, I discuss what is currently the best motivated 
and most completely developed physical model for the first moments of the universe: cosmological inflation 1, 2, 3 
Inflation naturally explains how the universe came to be so large, so old, and so flat, and provides a compellingly 
elegant and predictive mechanism for generating the primordial perturbations which gave rise to the rich structure 
we see in the universe today [71 |8j |9j [TOj [TTl [121 \lM ■ Inflation provides a link between the Outer Space of astrophysics 
and the Inner Space of particle physics, and gives us a window to physics at energy scales far beyond the reach of 
particle accelerators. Furthermore, inflation makes testable predictions, which have so far proven to be an excellent 
match to the data. 

The lectures are organized as follows: 

• Section [H provides an introduction and a brief overview of General Relativity. 

• Section [m discusses the Friedmann-Robertson- Walker spacetimc and the standard hot Big Bang picture of 
cosmology, including the Cosmic Microwave Background. 

• Section [TTTl explains unresolved issues in the standard cosmology, in particular the horizon and flatness problems. 

• Section llVI introduces inflation in scalar fleld theories. 

• Section [Vl discusses quantum fluctuations in inflation and the generation of cosmological perturbations. 

• Section [VI| discusses the observational predictions of inflation, and current constraints from Cosmic Microwave 
Background measurements. 

• Section [VIII discusses conclusions and the future outlook for inflationary physics. 

• Appendix \A\ describes in detail the generation of density perturbations during inflation. 



'Electronic address: |whkinney@buffal o.edu| 

^ Inflation in its current torm was introduced by Guth, but similar ideas had been discussed before [1][5]. A short history of the early 
development of inflation can be found in Ref. i6j. 
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The lectures are at an advanced undergraduate or beginning graduate student level. Most of the the lectures should 
be accessible with only a background in Special Relativity. A working knowledge of General Relativity and quantum 
field theory are helpful for Sections |IV| and |V] and for Appendix |X] Where possible, I reference review articles for 
further reading on related topics. For other reviews on inflation, see Refs. [51 [HI [T51 [T51 [T71 [T51 [TOl EU] . 



A. The Metric 



The fundamental object in General Relativity is the metric, which encodes the shape of the spacetime. A metric is 
a symmetric, bilinear form which defines distances on a manifold. For example, we can express Pythagoras' theorem 
in a Euclidean three-dimensional space, 

f ^x^ + y^ + z^, (1) 
as a matrix product over the identity matrix 6ij = diag (1, 1, 1), 

^^3^'^^- (2) 

Therefore the identity matrix 6ij can be identified as the metric for the Euclidean space: if we wish to describe a 
non- Euclidean manifold, we replace Sij with a more complicated matrix g^ , which in general can depend on the 
coordinates x'. For an arbitrary path through the space, we express distances on the manifold in differential form. 



df = Yg^jdx'dx'. (3) 



The distance along any path in the spacetime, or world line, is then given by integrating d£ along that path. A 
familiar example of a non-Euclidean space frequently used in physics is the Minkowski Space describing spacetime in 
Special Relativity. Distances along a world line in Minkowski Space are measured by the proper time, which is the 
time as measured by an observer traveling on that world line. The proper time s along a world line is given by the 
relation 

ds'^ = dt"^ ~ d-x^ 

= V^..dx''dx\ (4) 

/I, I/— 0,3 

where we take the speed of light c — 1. We express four-vectors as i = (t,x,y,z) — (x'^ , x^ , x'^ , x^) , and dx^ = 
dx^ + dy^ + dz^ is the Euclidean distance along a spatial interval. The metric rj^i, for Minkowski Space is given by 



1 \ 

-1 

-1 

-I I 



(5) 



Anything traveling the speed of light has velocity d\x.\ jdt — 1. Photons therefore always travel along world lines 
of zero proper time, ds^ — dt^ — dx.^ ~ 0, called null geodesies. Massive particles travel along world lines with real 
proper time, ds'^ > 0, called timelike geodesies. Causally disconnected regions of spacetime are separated by spacelike 
intervals, with ds^ < 0. The set of all null geodesies passing through a given point (or event) in spacetime is called 
the light cone (Fig. [ij The interior of the light cone, consisting of all null and timelike geodesies, defines the region 
of spacetime causally related to that event. 



B. General Relativity and the Einstein Field Equation 



The Minkowski metric rj^i, of Special Relativity describes a Euclidean spacetime which is static, empty, and infinite 
in space and time. The addition of gravity to the picture requires General Relativity, which describes gravitational 
fields as curvature in the spacetime. The fundamental object in General Relativity is the metric g^j^^^x), which 
describes the shape of the spacetime and in general depends on the spacetime coordinate x. As in Minkowski Space, 
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lengths in curved spacetime are measured by the proper time s, with the proper time along a world line determined 
by the metric 

ds'^= 9nu{x)dx''dx''. (6) 

^,1^—0,3 

As in Special Relativity, photons travel along null geodesies, with ds"^ = 0, and massive particles travel along timelike 
geodesies, with ds'^ > 0. However, unlike Special Relativity, null geodesies need not always be 45° lines defining light 
cones, but can be curved by gravity. 

In General Relativity, the distribution of mass/energy in the spacetime determines the shape of the metric, and the 
metric in turn determines the evolution of the mass /energy. Electromagnetism provides a convenient analogy: in elec- 
tromagnetism, the distribution of charges and currents determines the electromagnetic field, and the electromagnetic 
field in turn determines the evolution of the charges and currents. Given a current four- vector J^, Maxwell's Equations 
are a set of linear, first-order partial differential equations that allow us to calculate the resulting electromagnetic field 

A-TT 

d^F^"" = y d^F^"" = — J*^. (7) 

c 

L/=0,3 

Here we have explicitly included the speed of light c to highlight its role as an electromagnetic coupling constant. 

We also adopt the typical summation convention for relativity: repeated indices are implicitly summed over. In 
General Relativity, we describe the distribution of mass/energy in a covariant way by specifying a symmetric rank-2 
stress-energy tensor T^j^y, which acts as a source for the gravitational field similar to the way the current four- vector 
sources electromagnetism. The analog of Maxwell's Equations is the Einstein Field Equation, which can be written 
in the deceptively simple form 



(8) 
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where the coupling constant is Newton's gravitational constant G. The tensor G^ijy, called the Einstein Tensor, is a 
symmetric 4x4 tensor consisting of the metric and its first and second derivatives. The Einstein Field Equation 
therefore represents a set of ten coupled, nonlinear, second-order partial differential equations of ten free functions, 
which are the elements of the metric tensor g^^. However, only six of these equations are actually indpendent, leaving 
four degrees of freedom. The physics of gravity is independent of coordinate system, and the additional degrees of 
freedom correspond to a choice of a coordinate system, or gauge on the four-dimensional space. Gravity is much more 
complicated than electromagnetism! As with any intractably complicated problem, we simplify the job by introducing 
a symmetry. In General Relativity there are a number of symmetries which allow either exact or perturbative solution 
to the Einstein Field Equations: 

• Vacuum: T^^ = 0. If we evaluate the Einstein Field Equations for small perturbations about an empty Minkowski 
Space, we find that they reduce at lowest order to a wave equation, and therefore General Relativity predicts 
the existence of gravity waves. 

• Spherical Symmetry. If we assume a spherically symmetric spacetime (also empty of matter, T^^ = 0) the 
Einstein Field Equation can be solved exactly, resulting in the Schwarzschild solution for black holes. 

• Homogeneity and Isotropy. If we assume that the stress-energy is distributed in a fashion which is homogeneous 
and isotropic, this is called a Friedmann- Roberts on- Walker (FRW) space, and is the case of interest for cosmology. 
Since the homogeneity and isotropy remove all spatial dependence, the Einstein Field Equations reduce from a 
set of partial differential equations to a set of nonlinear ordinary differential equations in time. For particular 
types of homogeneous, isotropic matter, these equations can be solved exactly, and perturbations about those 
exact solutions can be handled self-consistently. 

Continuing the analogy with electromagnetism, the equivalent of charge conservation, 

d^J'^ = ^ + V-i=0, (9) 
in General Relativity is stress-energy conservation 

= 0, (10) 

where I?^ represents a covariant derivative, which is a generalization of the partial derivative to a curved manifold. 
We will also denote covariant derivatives with a semicolon, for example T'^''.^, — 0. Likewise, simple partial derivatives 
are denoted with a comma, df /dx^ = 9^/ = /_^. As in the case of electromagnetism, where the charge conservation 
equation is not independent, but is instead a consequence of Maxwell's Equations, stress-energy conservation in General 
Relativity is a consequence of the Einstein Field Equations and does not independently constrain the solutions. In 
the next section, we discuss FRW spaces and their application to cosmology in more detail. 



II. FRIEDMANN-ROBERTSON- WALKER SPACETIMES 



A. The Friedmann Equation 



A homogeneous space is one which is translationally invariant, or the same at every point. An isotropic space is one 
which is rotationally invariant, or the same in every direction. The two are not the same: a space which is everywhere 
isotropic is necessarily homogeneous, but a space which is homogeneous is not necessarily isotropic. (Consider, for 
example a space with a uniform electric field: it is translationally invariant but not rotationally invariant.) It is possible 
to show |21j that the most general metric consistent with homogeneity and isotropy is obtained by multiplying a static 
spatial geometry with a time-dependent scale factor a{t): 



ds^ 



dt^ 
dt^ 



{t) dx^ 



dr^ 



1 — kr^ 



r'^dn' 



(11) 



where we have expressed the spatial line element in terms of spherical coordinates r,9,(f), and the solid angle is given 
by the usual dO^ = s'm OdOdcf). The constant k defines the curvature of the spacetime, with fc = corresponding 
to flat (Euclidean) spatial sections, and fc = ±1 corresponding to positive and negative curvatures, respectively. A 



5 



spacetime of this general form is called a Friedmann- Roberts on- Walker (FRW) spacetime. Likewise, the most general 
homogeneous, isotropic stress-energy is diagonal, with all of its spatial components identical, 



/Pit) 



V 




(12) 



where we identify the energy density p and the pressure p from the continuity equation arising from stress-energy 
conservation, 

T^% = P + 3(^iyp + p)=0. (13) 
The Einstein field equations then reduce to a set of two coupled, non-linear ordinary differential equations, 



2 

a ' 



k _ Sir 

a\ 47r , , , , 

-Tr^(p + 3p). (14) 



The first is called the Friedmann Equation, and the second is called the Raychaudhuri Equation. Note that the 
equations for the evolution of the scale factor depend not only on the energy density p, but also the pressure p: 



pressure gravitates! The continuity equation (13 1 is not independent of the Einstein Field Equations (14 1, but can 
be derived directly from the Friedmann and Raychaudhuri Equations. The expansion rate a/ a is called the Hubble 
parameter H: 

H^^, (15) 
a 

and has units of inverse time. A positive Hubble parameter H > corresponds to an expanding universe, and a 
negative Hubble parameter H < corresponds to a collapsing universe. (Since our actual universe is expanding, we 
will specialize to that case.) Minkowski Space can be recovered by assuming a flat geometry k — 0, and no expansion, 
d = 0. The Hubble parameter sets the fundamental scale of the spacetime, i.e. a characteristic time is the Hubble 
time t ~ H~^, and likewise the Hubble length is c? ~ H~^. We will see later that the Hubble time sets the scale for 
the age of the universe, and the Hubble length sets the scale for the size of the observable universe. 

The coordinate system (t, x) is called a comoving coordinate system, because observers with constant comoving 
coordinates are at rest relative to the expansion, i.e. two observers with constant separation in comoving coordinates 
Ax have a physical, or proper, separation which increases in proportion to the scale factor 

Axpi-op = a (t) Axcom- (16) 

An important kinematic effect of cosmological expansion is the phenomenon of cosmological redshift: we will see later 
that solutions to the wave equation in an FRW space have constant wavelength in comoving coordinates, so that the 
proper wavelength of (for example) a photon increases in time in proportion to the scale factor 

Acxa(t). (17) 

For a photon emitted at time tom and detected at time to, the redshift z is defined by: 

(l + ^)-T^ = ^. (18) 

(Here we introduce the convention used frequently in cosmology that a subscript refers to the current time, not 
an initial time.) Note that the cosmological redshift is not a Doppler shift caused by the relative velocity of the 
source and detector, but is an expansion effect: the wavelength of a photon traveling through the spacetime increases 
because the underlying spacetime is expanding. Another way to look at this is that a photon traveling through an 
FRW spacetime loses momentum with time, 

p= hiy (X a~^{t). (19) 
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By the equivalence principle, this momentum loss must apply to massive particles as well as photons: any particle 
moving in an expanding FRW spacetime will lose momentum as p cx a~^. For massless particles like photons, this is 
manifest as a redshift in the wavelength, but it means that a massive particle will asymptotically come to rest relative 
to the comoving coordinate system. Thus, comoving coordinates represent a preferred reference frame reminiscent of 
Aristotelian physics: any free body with a "peculiar" velocity relative to the comoving frame will eventually come to 
rest in that frame. 

There are three possibilities for the curvature of the universe: flat (k = Q), positively curved (fc = +1), or negatively 
curved {k = —1). The current value of the Hubble parameter is (from the Hubble Space Telescope Key Project [22), 

i7o = 72 ± 8 km/s/Mpc. (20) 



Therefore, we can see from the Friedmann Equation ( 14 1 that, given the expansion rate the curvature is determined 
by the density: 

''-'[■^r"')- '''' 

Note that only the sign of k is physically important, since any rescaling of k is equivalent to a rescaling of the scale 
factor a. We define the critical density as the density for which fc = 0, corresponding to a geometrically flat universe, 

^ fc = o. (22) 

For p > pc the universe is positively curved and closed, with finite volume, and for p < pc, the universe is negatively 
curved and open, with infinite volume. We express the ratio of the actual density p to the critical density pc as the 
parameter fi: 

(Do not confuse the density parameter fl with the solid angle dfl in Eq. 11 ) Table [l] summarizes the relation between 
density, curvature, and geometry. The density parameter il is not in general constant in time, and we can re-write 



TABLE I: Cosmological density and curvature 



density 


curvature 


geometry 


n = 1 


k = 


flat 


O > 1 


k = 1 


closed 


n < 1 


k = 


open 



the Friedmann Equation as 

^i*) = ^ + 7^- (24) 
(aH) 

Since the Hubble parameter is proportional to the inverse time H cx t~^, we see that the time-dependence of fl is 
determined by the time dependence of the scale factor a (t). In the next section, we tackle the problem of solving for 
a{t). 



B. Solving the Friedmann Equation 

In the previous section, we considered the form and kinematics of FRW spaces, but not the dynamics, that is, how 
does the stress-energy of the universe determine the expansion history? The answer to this question depends on what 
kind of matter dominates the cosmological stress-energy. In this section, we consider three basic types of cosmological 
stress-energy: matter, radiation, and vacuum. 

The simplest kind of cosmological stress-energy is generically referred to as matter. Imagine a comoving box with 
sides of length L. By comoving box, we mean a box whose corners are at rest in a comoving coordinate system, 
and whose proper dimension is therefore increasing proportional to the scale factor, Lprop oc a. That is, the box is 
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FIG. 3: A comoving box full of radiation. The number density of photons scales inversely with the volume of the box, but the 
photons also increase in wavelength. 

growing with the expansion of the universe. Now imagine the box filled with N particles of mass m, also at rest in the 
comoving reference frame (Fig. [2]). In units where c = 1, the relativistic energy density of such a system of particles 
is given by 




FIG. 4: A comoving box full of vacuum. The energy density of vacuum does not scale at all! 
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where V is the proper volume of the box, V = ipiop a^. Since neither M nor N change with expansion, we have 
immediately that 

Pm = 7^ oc a , (26) 

where L is the comoving size of the box. So the proper energy density of massive particles at rest in a comoving 
volume evolves as the inverse cube of the scale factor. Now imagine the same box filled with N photons with frequency 
v (Fig. [3|. The energy per photon is hv, so that the energy density in the box is then 

Pi = —- (27) 

As in the case of massive particles, the number density of photons in the box redshifts inversely with the pr ope r 
volume of the box n — N/V oc a^^. But each photon also loses energy through cosmological redshift, i' cx a^^ (19 1, 



so that the total energy density in photons or other massless degrees of freedom, which we generically refer to as 
radiation, redshifts as 

cx a"^. (28) 

Note also that cosmological redshift immediately gives us a rule for the behavior of a black-body spectrum of radiation 
with temperature T. Since all photons redshift at exactly the same rate, a system with starts out as a black-body 
stays a black-body, with a temperature that decreases with expansion, 

T^cxa"^ (29) 

The third type of stress-energy which is important in cosmology dates back to Einstein's introduction of a "cosmo- 
logical constant" to his field equations. If we take the stress-energy T^^^ and add a term proportional to the metric. 



the identity Di^g^^" — means the stress-energy conservation equation ( 10 ) is unchanged 



D^T'"' ^ D^iTf"" + Ag'"')^0. (30) 

In our analogy with electromagnetism, this is like adding a constant to the electromagnetic potential, V'{x) — V{x)+A. 
The constant A does not affect local dynamics in any way, but it does affect the cosmology. From Eq. ( [T2| ), stress- 
energy of the form T^'^ — Kg^'^ corresponds to an equation of state 

PA = -/OA- (31) 



The continuity equation (13 1 then reduces to 



P+3(^^j (P + P) = P=0, (32) 

so that vacuum has a constant energy density, pA ~ const. A cosmological constant is also frequently referred to 
as vacuum energy, since it is as if we are assigning an energy density to empty space. With this interpretation, a 
comoving box full of vacuum contains a total amount of energy which grows with the expansion of the universe (Fig. 
|4]). This highlights the curious property of General Relativity that, while energy is conserved in a local sense, it is 
not conserved globally. We are creating energy out of nothing! 

It is straightforward to solve the Einstein Field Equations for the three basic types of stress-energy. Consider first 
a matter-dominated universe. We can write the time derivative of the energy density as: 

p,n cx a^^ ^ pm = -3 f ^ j p. (33) 



From the continuity equation (13 1, we have 



P+3(^j(p + .)^3(^j,^0. (34) 
We then have that the pressure of matter vanishes, Pm = 0. The matter-dominated Friedmann Equation becomes 



k 



a } 3TOpj 



p cx a"^. (35) 
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FIG. 5: Schematic diagram of how the three types of stress-energy scale with redshift 1 + z <x a: at early time, radiation 
dominates, followed by matter, and finally the universe is dominated by vacuum energy. 



In the case of a flat universe, k — 0, the solution is especially simple 

• \ 2 

a 
a 



-) ^a{t)<xt^^\ (36) 



Similarly, for a radiation dominated universe, the continuity equation implies that 

oc a^^ ^ — p-y/3. (37) 



Again assuming a flat geometry. 



• \ 2 

a 
a 



cx a^^ ^ a{t) (X i^/^ (38) 



Finally, solving the the Friedmann Equation for the vacuum case gives 

2 



X PA = const. ^ a (t) (X e , (39) 
so that the universe expands exponentially quickly, with a time constant given by the Hubble parameter 

i/=^— PA = const. (40) 

Such a spacetime is called de Sitter space. 

Note in particular that the energy density in radiation redshifts away more quickly than the energy density in 
matter, and vacuum energy does not redshift at all, so that a universe with a mix of radiation, matter and vacuum 
will be radiation-dominated at early times, matter-dominated at later times, and eventually vacuum-dominated (Fig. 
[5]). Note also that for either matter- or radiation-domination, the universe is singular as t — s- 0: the universe has finite 
age! Since the scale factor vanishes at i = 0, and the density scales as an inverse power of a, the initial singularity 
consists of infinite density. Likewise, since temperature also scales inversely with a, the initial singularity is also a 
point of infinite temperature. We therefore arrive at the standard hot Big Bang picture of the universe: a cosmological 
singularity at finite time in the past, followed by a hot, radiation-dominated expansion, during which the universe 
gradually cools as T oc and the radiation dilutes, followed by a period of matter-dominated expansion during 
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which galaxies and stars and planets form. Finally, if the vacuum energy is nonzero, it will inevitably dominate, and 
the universe will enter a state of exponential expansion. Current evidence indicates that the real universe made a 
transition from matter-domination to vacuum-domination at a redshift of around 2; = 1, or about a billion years ago, 
so that the densities of the three types of matter today are of order 





-0.7, 




~ 0.3, 


fly 


~ 10" 



In the next section, we discuss one important prediction of the hot Big Bang: the presence of a background of relic 
photons from the early universe, called the Cosmic Microwave Background. 



C. The Hot Big Bang and the Cosmic Microwave Background 

The basic picture of an expanding, cooling universe leads to a number of startling predictions: the formation of nuclei 
and the resulting primordial abundances of elements, and the later formation of neutral atoms and the consequent 
presence of a cosmic background of photons, the Cosmic Microwave Background (CMB) [23j EH EH |26l [21] . A rough 
history of the universe can be given as a time line of increasing time and decreasing temperature [5S]: 

• T = 00, t = 0: Big Bang. 

• T ^ 10^^ K, t ^ 10"^^ sec: Primordial soup of fundamental particles. 

• T ^ 10^'^ K, t ^ 10"^ sec: Protons and neutrons form. 

• T ^ lO^'^ K, t ^ ?! min: Nucleosynthesis: nuclei form. 

• T - 3000 K, t - 300, 000 years: Atoms form. 

• T ^ 10 i^, t ~ 10^ years: Galaxies form. 

• T - 3 if , t - 10^° years: Today. 

The epoch at which atoms form, when the universe was at an age of 300,000 years and at a temperature of around 
3000 K is oxy moronic ally referred to as "recombination" , despite the fact that electrons and nuclei had never before 
"combined" into atoms. The physics is simple: at a temperature of greater than about 3000 K, the universe consisted 
of an ionized plasma of mostly protons, electrons, and photons, with a few helium nuclei and a tiny trace of lithium. 
The important characteristic of this plasma is that it was opaque, or, more precisely, the mean free path of a photon 
was a great deal smaller than the Hubble length. As the universe cooled and expanded, the plasma "recombined" 
into neutral atoms, first the helium, then a little later the hydrogen. 

If we consider hydrogen alone, the process of recombination can be described by the Saha equation for the equilibrium 
ionization fraction of the hydrogen [25] : 

1-Xc 4\/2C(3) ( tV'^ /13.6eV\ 

- —ni — ] exp — ^ . (42) 



Here rric is the electron mass and 13.6 eV is the ionization energy of hydrogen. The physically important parameter 
affecting recombination is the density of protons and electrons compared to photons. This is determined by the haryon 
asymmetry,^ which is described as the ratio of baryons to photons: 

= ^^^^^^ - 2.68 x 10-« in^h^) . (43) 

Here fib is the baryon density and h is the Hubble constant in units of 100 km/s/Mpc, 

h = Ho/{lQQ km/s/Mpc). (44) 



^ If there were no excess of baryons over antibaryons, there would be no protons and electrons to recombine, and the universe would be 
just a gas of photons and neutrinos! 
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FIG. 6: Schematic diagram of recombination. At early time, the temperature of the universe is above the ionization energy of 
hydrogen and helium, so that the universe is full of an ionized plasma, and the mean free path for photons is short compared 
to the Hubble length. At late time, the temperature drops and the nuclei capture the electrons and form neutral atoms. Once 
this happens, the universe becomes transparent to photons, which free stream from the surface of last scattering. 



The most recent result from the WMAP satellite gives Uy^h^ = 0.02273±0. 00062 [29 . Recombination happens quickly 
(i.e., in much less than a Hubble time t ~ H~^), but it is not instantaneous. The universe goes from a completely 
ionized state to a neutral state over a range of redshifts Az ~ 200. If we define recombination as an ionization fraction 
Xg = 0.1, we have that the temperature at recombination = 0.3 eV. 

What happens to the photons after recombination? Once the gas in the universe is in a neutral state, the mean free 
path for a photon becomes much larger than the Hubble length. The universe is then full of a background of freely 
propagating photons with a blackbody distribution of frequencies. At the time of recombination, the background 
radiation has a temperature of T = Tr = 3000 K, and as the universe expands the photons redshift, so that the 
temperature of the photons drops with the increase of the scale factor, T cx a{t)~^. We can detect these photons 
today. Looking at the sky, this background of photons comes to us evenly from all directions, with an observed 
temperature of Tq ~ 2.73 K. This allows us to determine the redshift of recombination. 



1 + ^R 



a{tn) 



Tr 
To 



1100. 



(45) 



This is the cosmic microwave background. Since by looking at higher and higher redshift objects, we are looking 
further and further back in time, we can view the observation of CMB photons as imaging a uniform "surface of last 
scattering" at a redshift of 1100 (Fig. [?]). 

To the extent that recombination happens at the same time and in the same way everywhere, the CMB will be 
of precisely uniform temperature. While the observed CMB is highly isotropic, it is not perfectly so. The largest 
contribution to the anisotropy of the CMB as seen from earth is simply Doppler shift due to the earth's motion through 
space. (Put more technically, the motion is the earth's motion relative to a comoving cosmological reference frame.) 
CMB photons are slightly blueshifted in the direction of our motion and slightly redshifted opposite the direction of 
our motion. This blueshift/redshift shifts the temperature of the CMB so the effect has the characteristic form of a 
"dipole" temperature anisotropy (Fig. [s]). The dipole anisotropy, however, is a ZocaZ phenomenon. Any intrinsic, or 
primordial, anisotropy of the CMB is potentially of much greater cosmological interest. To describe the anisotropy of 
the CMB, we remember that the surface of last scattering appears to us as a spherical surface at a redshift of 1100. 
Therefore the natural parameters to use to describe the anisotropy of the CMB sky is as an expansion in spherical 




FIG. 7: Cartoon of the last scattering surface. From earth, we see blackbody radiation emitted uniformly from all directions, 
forming a "sphere" at redshift z = 1100. 



harmonics 

AT °° ^ 

-Y = Y.Y. ai„,Y,m{e,4>). (46) 

1=1 m=-t 

If we assume isotropy, there is no preferred direction in the universe, and we expect the physics to be independent of 
the index m. We can then define 



^El«^™l'- (47) 



The i = 1 contribution is just the dipole anisotropy. 



10"^ (48) 



AT 

The dipole was first measured in the 1970 's by several groups [BO] , 131 ] [32] . It was not until more than a decade after 
the discovery of the dipole anisotropy that the first observation was made of anisotropy for ^ > 2, by the differential 
microwave radiometer aboard the Cosmic Background Explorer (COBE) satellite [33], launched in in 1990. COBE 
observed that the anisotropy at the quadrupole and higher t was two orders of magnitude smaller than the dipole: 



Ar\ 



~ 10"^ (49) 



Fig. [8] shows the dipole and higher-order CMB anisotropy as measured by COBE. This anisotropy represents 
intrinsic fluctuations in the CMB itself, due to the presence of tiny primordial density fluctuations in the cosmological 
matter present at the time of recombination. These density fluctuations are of great physical interest, since these 
are the fluctuations which later collapsed to form all of the structure in the universe, from superclusters to planets 
to graduate students. While the physics of recombination in the homogeneous case is quite simple, the presence of 
inhomogeneities in the universe makes the situation much more complicated. I describe some of the major effects in 
a qualitative way here, and refer the reader to the literature for a more detailed technical explanation of the relevant 
physics [231 im [55J |57]. In these lectures, I primarily focus on the current status of the CMB as a probe of 
inflation, but there is much more to the story. 
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FIG. 8: The COBE measurement of the CMB anisotropy [33] • The top oval is a map of the sky showing the dipole anisotropy 
AT/T ~ 10"^. The bottom oval is a similar map with the dipole contribution and emission from our own galaxy subtracted, 
showing the anisotropy for £ > 1, AT/T ~ 10~^. (Figure courtesy of the COBE Science Working Group.) 




FIG. 9: The WMAP measurement of the CMB anisotropy [34]. (Figure courtesy of the WMAP Science Working Group.) 
WMAP measured the anisotropy with much higher sensitivity and resolution than COBE. 



14 



The simplest contribution to the CMB anisotropy from density fluctuations is just a gravitational redshift, known 
as the Sachs-Wolfe effect [35] , A photon coming from a region which is slightly denser than the average will have 
a shghtly larger redshift due to the deeper gravitational well at the surface of last scattering. Conversely, a photon 
coming from an underdense region will have a slightly smaller redshift. Thus we can calculate the CMB temperature 
anisotropy due to the slightly varying Newtonian potential $ from density fluctuations at the surface of last scattering: 

ST 1 

^ = 3 [<f>cm - 1>obs] , (50) 

where <i>oin is the potential at the point the photon was emitted on the surface of last scattering, and <i>obs is the 
potential at the point of observation, which can be treated as a constant. The factor 1/3 is a General Relativistic 
correction. This simple kinematic contribution to the CMB anisotropy is dominant on large angular scales, corre- 
sponding to multipoles £ < 100. However, the amount of information we can gain from these multipoles is limited 
by an intrinsic source of error called cosmic variance. Cosmic variance is a result of the statistical nature of the 
primordial power spectra: since we have only one universe to measure, we have only one realization of the random 
fleld of density perturbations, and therefore there is an inescapable I/Vn uncertainty in our ability to reconstruct 
the primordial power spectrum, where N is the number of independent wave modes which will flt inside the horizon 
of the universe! On very large angular scales, this problem becomes acute, and we can write the cosmic variance error 
on any given Ci as 

^ - ' (51) 

which comes from the fact that any Ci is represented by 2£+l independent amplitudes a^^. Even a perfect observation 
of the CMB can only approximately measure the true power spectrum — the errors in the WMAP data, for example. 



are dominated by cosmic variance out to ^ ~ 400 (Fig. 10) 



For fluctuation modes on smaller angular scales, more complicated physics comes into play. The dominant process 
that occurs on short wavelengths is acoustic oscillations in the baryon/photon plasma. The idea is simple: matter 
tends to collapse due to gravity onto regions where the density is higher than average, so the baryons "fall" into 
over dense regions. However, since the baryons and the photons are still strongly coupled, the photons tend to 
resist this collapse and push the baryons outward. The result is "ringing" , or oscillatory modes of compression and 
rarefaction in the gas due to density fluctuations. The gas heats as it compresses and cools as it expands, which 
creates fluctuations in the temperature of the CMB. This manifests itself in the Ci spectrum as a series of peaks and 
valleys (Fig. [lO]). The specific shape and location of the acoustic peaks is created by complicated but well-understood 
physics, involving a large number of cosmological parameters. The presence of acoustic peaks in the CMB was first 
suggested by Sakharov , and later calculated by Sunyaev and Zel'dovich |371 |3H] and Peebles and Yu [31] . The 
complete linear theory of CMB fluctuations was worked out by Ma and Bertschinger in 1995 |1D] . The shape of the 
CMB multipole spectrum depends, for example, on the baryon density fib, the Hubble constant Hq, the densities 
of matter and cosmological constant fl^, the amplitude of primordial gravitational waves, and the redshift 
at which the first generation of stars ionized the intergalactic medium. This makes interpretation of the spectrum 
something of a complex undertaking, but it also makes it a sensitive probe of cosmological models. 

In addition to anisotropy in the temperature of the CMB, the photons coming from the surface of last scattering 
are expected to be weakly polarized due to the presence of perturbations [411 142] . This polarization is much less 
well measured than the temperature anisotropy, but it has been detected by WMAP and by a number of ground- 
and balloon-based measurements [43, 44, 45) UHl H?]. Measurement of polarization promises to greatly increase the 
amount of information it is possible to extract from the CMB. Of particular interest is the odd-parity, or B-mode 
component of the polarization, the only primordial source of which is gravitational waves, and thus provides a clean 
signal for detection of these perturbations. The B-mode has yet to be detected by any measurement. 



III. THE FLATNESS AND HORIZON PROBLEMS 



We have so far considered two types of cosmological mass-energy - matter and radiation - and solved the Friedmann 
Equation for the simple case of a fiat universe. What about the more general case? In this section, we consider non-fiat 
universes with general contents. We introduce two related questions which are not explained by the standard Big 
Bang cosmology: why is the universe so close to fiat today, and why is it so large? 

We can describe a general homogeneous, isotropic mass-energy by its equation of state 



p = wp, 



(52) 
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FIG. 10: The Ci spectrum for the CMB as measured by WMAP, showing the peaks characteristic of acoustic oscillations. The 
gray shaded region represents the uncertainty due to cosmic variance. (Figure courtesy of the WMAP Science Working Group.) 



so that pressureless matter corresponds to w — Q, and radiation corresponds to w = 1/3. We will consider only the 
case of constant equation of state, w = const. From the continuity equation, we have 



p+3(l + u;)-p = 0, 



with solution 



p (X a 



The Friedmann Equation for a flat universe is then 



-3{l+w) 



-3(l+«;) 



SO that the scale factor increases as a power-law in time, 

a(t)oci2/3(i+-). 



(53) 



(54) 



(55) 



(56) 



What about the evolution of a non-flat universe? Analytic solutions for a{t) in the k ^ case can be found in 
cosmology textbooks. For our purposes, it is sufficient to consider the time-dependence of the density parameter f2. 
From Eqs. (14 23 24) it is not too difficult to show that the density parameter evolves with the scale factor a as: 



din a 



{i + 3w)n{n-i). 



(57) 



Proof is left as an exercise for the reader. Note that a flat universe, = 1 remains flat at all times, but in a non-flat 
universe, the density parameter is a time-dependent quantity, with the evolution determined by the equation of 
state parameter w. For matter (w — 0) or radiation (w — 1/3), the prefactor in Eq. (57 1 is positive. 



1 + 3w > 0, 

which means a flat universe is an unstable fixed point: 

d\n- 1 



din a 



> 0, (1-1- 3w) > 0. 



(58) 



(59) 
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Any deviation from a flat geometry is amplified by the subsequent cosmological expansion, so a nearly flat universe 
today is a highly fine-tuned situation. The WMAP5 CMB measurement tells us the universe is flat to within a few 
percent, jfio ^ 1| < 0.02 [551 HH]- If we are very conservative and take a limit on the density today as fip = 1 ± 0.05, 
that means that at recombination, when the CMB was emitted, firoc = 1 ± 0.0004, and at the time of primordial 
nucleosynthesis, finuc = 1 ± 10^^^. Why did the universe start out so incredibly close to flat? The standard Big Bang 
cosmology provides no answer to this question, which we call the flatness problem. 

There is a second, related problem with the standard Big Bang picture, arising from the finite age of the universe. 
Because the universe has a finite age, photons can only have traveled a finite distance in the time since the Big Bang. 
Therefore, the universe has a horizon: the further out in space we look, the further back in time we see. If we could 
look far enough out in any direction, past the surface of last scattering, we would be able to see the Big Bang itself, 
and beyond that we can see no further. Every observer in an FRW spacetime sees herself at the center of a spherical 
horizon which defines her observable universe. To calculate the size of our horizon, we use the fact that photons travel 
on paths of zero proper length: 

ds^ = dt^ - {t) \diif = 0, (60) 
so that the comoving distance |(ix| traversed by a photon in time dt is 

|rfx| = (61) 
Therefore, the size of the cosmological horizon at time t after the Big Bang is 

Jo ait') 

To convert comoving length to proper length, we just multiply by a (t), so that the proper horizon size is 

dr''{t) = a{t)d^^"^t). (63) 

Normalizing a (to) — 1, the horizon size of a 14-billion year-old flat, matter-dominated universe is dn = Sto ^ 13 Gpc. 

To see why the presence of a horizon is a problem for the standard Big Bang, we examine the causal structure of 
an FRW universe. Take the FRW metric 

ds^ ^ dt^ - (t) \dyif , (64) 
and re-write it in terms of a redefined clock, the conformal time r: 

ds^ = (t) \dT^ - \dx\^] . (65) 



Conformal time is a "clock" which slows down with the expansion of the universe, 



dT=-^, (66) 
a(t) 



so that the comoving horizon size is just the age of the universe in conformal time 

/■* dt' r 

dn{t)= / / dr' = r. (67) 

The conformal metric is useful because the expansion of the spacetime is factored into a static metric multiplied by a 
time-dependent conformal factor a (r), so that photon geodesies are simply described by d |x| = dr. In a diagram of r 
versus |x|, photons travel on 45° angles. (Note that this is true even for curved spacetimes!) We can draw light cones 
and infer causal relationships with the expansion factored out, in a manner identical to the usual case of Minkowski 
Space. 

There is one major difference between FRW and Minkowski: an FRW spacetime has a finite age. Therefore, unlike 
the case of Minkowski Space, which has an infinite past, an FRW spacetime is "chopped off" at some finite past time 



T — (Fig 11 ). The initial singularity is a surface of constant conformal time, and it is easy to see from Eq. (67 1 that 
our horizon size is the width of our past light cone projected on the surface defined by the initial singularity. This is 
a very different picture from the notion many people (even scientists) have of the Big Bang, which is something akin 
to an explosion, with the universe initially a cosmic "egg" of zero size. On the contrary, in the case of a flat or open 
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FIG. 11: A conformal diagram of a Friedmann-Robertson- Walker space. The FRW space is causally identical to Minkowski 
Space, except that it is not past-infinite, so that past light cones are "cut off" at the Big Bang, which is a spatially infinite 
surface at time t = 0. 



universe, the universe is spatially infinite an infinitesimal amount of time after the initial singularity: the Big Bang 
happens everywhere at once in an infinite space! Our observable universe is finite because we can only see a small 
patch of the much larger cosmos.^ Closed universes are spatially finite, but are still much larger in extent than our 
observable patch. The key point is that two events on the conformal spacetime diagram are causally connected only 
if they share a causal past: that is, if their past light cones overlap. 



Consider two points on the CMB sky 180° degrees apart (Fig. 12 1. Their past light cones do not overlap, and the 
two points are causally disconnected. Those two points on the surface of last scattering occupy completely separate, 
disconnected observable universes. How did these points reach the observed thermal equilibrium to a few parts in 
10^ if they never shared a causal past? This apparent paradox is called the horizon problem: the universe somehow 
reached nearly perfect equilibrium on scales much larger than the size of any local horizon. From the Friedmann 
Equation, it is easy to show that the horizon problem and the fiatness problem are related: consider a comoving 
length scale A. It is is easy to show that for w = const., the ratio of A to the horizon size is related to the curvature 
by a conservation law 



const. 



Proof is left as an exercise for the reader. Therefore, for a universe evolving away from flatness, 

d\n-i\ 



din a 



>0, 



(68) 



(69) 



Of course, this is an idealization, and the actual universe could well have a nontrivial global topology, even if it is locally flat, as long 
as the scale of the overall manifold is much larger than our horizon size [491 150) . 
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FIG. 12: A conformal diagram of the Cosmic Microwave Background. Two points on opposite sides of the sky are causally 
separate, since their past hght cones do not intersect. 



the horizon size gets bigger in comoving units 



A 



d In a V dfj 



< 0. (70) 



That is, more and more space "falls into" the horizon, or becomes causally connected, at late times. 

What would be required to have a universe which evolves toward flatness, rather than away from it? From Eq. 



( 57 1 , we see that having 1 + 3w negative will do the trick, 

< 0, (1 + S^i;) < 0. (71) 

din a 

Therefore, if the energy density of the universe is dominated not by matter or radiation, but by something with 
sufficiently negative pressure, p < — p/3, a curved universe will become flatter with time. From the Raychaudhuri 



Equation (14 1, we see that the case of p < — p/3 is exactly equivalent to an accelerating expansion: 

- oc -(l + Sw) > 0, (l + 3w) < 0. (72) 
a 

If the expansion of the universe is slowing down, as is the case for matter- or radiation-domination, the curvature 



evolves away from flatness. But if the expansion is speeding up, the universe gets flatter. From Eq. (68 1, we see that 
this negative pressure solution also solves the horizon problem, since accelerating expansion means that the horizon 
size is shrinking in comoving units: 



d /A 
c? In a \ fffj 



> 0, (1 + 3w) < 0. (73) 



When the expansion accelerates, distances initially smaller than the horizon size are "redshifted" to scales larger than 
the horizon at late times. Accelerating cosmological expansion is called inflation. 
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The simplest example of an accelerating expansion from a negative pressure fluid is the case of vacuum energy we 



considered in Section ( II B I , for which the scale factor increases exponentially, 



a (X e 



Ht (74) 



For such expansion, the universe is driven exponentially toward a flat geometry, 



^^""=2(1-17). (75) 



din a 

We can see that the horizon problem is also solved by looking at the conformal time: 

dt 



dr^^^ e-"'dt, (76) 
a{t) 



so that 



The conformal time during the inflationary period is negative, tending toward zero at late time. Therefore, if we have 
a period of inflationary expansion prior to the early epoch of radiation-dominated expansion, inflation takes place in 
negative conformal time, and conformal time t — Q represents not the initial singularity but the transition from the 
inflationary expansion to radiation domination. The initial singularity is pushed back into negative conformal time, 
and can be pushed arbitrarily far depending on the duration of inflation. Figure [13] shows the causal structure of an 
inflationary spacetime. The past light cones of two points on the CMB sky do not intersect at t — 0, but inflation 
provides a "sea" of negative conformal time, which allows those points to share a causal past. In this way, inflation 
solves the horizon problem. 

In more realistic models of inflation in the early universe, the energy density is approximately, but not exactly, 
constant, and the expansion is approximately, but not exactly, exponential. In such quasi-de Sitter spaces, the quali- 
tative picture above still holds, and inflation provides a clean and compelling explanation for the peculiar boundary 
conditions for our universe. In the next section, we discuss how to construct more detailed models of inflation in fleld 
theory. 



IV. INFLATION FROM SCALAR FIELDS 



The example of de Sitter evolution we considered in Section |III| gives a good qualitative picture of how inflation, 
or accelerated expansion, solves the horizon and flatness problems of the standard Big Bang cosmology. However, 
this leaves open the question: what physics is responsible for the accelerated expansion at early times? It cannot 
be Einstein's cosmological constant, simply because a universe dominated by vacuum energy stays dominated by 
vacuum energy for the infinite future, since in a de Sitter background matter (p oc a~^) and radiation (p cx a~^) are 
diluted exponentially quickly. Therefore, we will never reach a radiation-dominated phase, and we will never see a 
hot Big Bang. In order to transition from an inflating phase to a thermal equilibrium, radiation-dominated phase, 
the vacuum-like energy during inflation must be time-dependent. We model this dynamics with a scalar field (/>, for 
which we assume the following action: 



where g = Det {g^i,) is the determinant of the metric and the Lagrangian for the field is 

^0 - \g^''d^<j,d,(t> ~V{(t>). 



(78) 



(79) 



a classical field theory to curved spacetime: 

'^'Minkowski 



Comparing the action ( 78 1 and the Lagrangian ( 79 1 with their Minkowski counterparts illustrates how we generalize 

d^x 



(80) 



The metric appears in two places in the curved-spacetime action: First, it appears in the measure of volume in 
the four-space, d'^x, where the determinant of the metric takes the role of the Jacobian for arbitrary coordinate 
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FIG. 13: A conformal diagram of light cones in an inflationary universe. Inflation ends in reheating at conformal time t — 0, 
which is the onset of the radiation-dominated expansion of the hot Big Bang. However, inflation provides a "sea" of negative 
conformal time, which allows the past light cones of events at the last scattering surface to overlap. 



transformations, 



Second, the metric appears in the kinetic term for the scalar field, where we replace the 



Minkowski metric 77'"^ with the general metric g'"'. 

The action (78 1 is not the most general assumption we could make, as we can see by writing the foil action including 
gravity. 



IGtt 



(81) 



Here R is the Ricci Scalar, composed of the metric and its derivatives. Variation of the first term in the action 
results in the Einstein Field Equation ([8|. Such a minimally coupled theory assumes that there is no direct coupling 
between the field and the metric, which would be represented in a more general action by terms which mix R and 
(j). In practice, many such non-minimally coupled theories can be transformed to a minimally coupled form by a field 
redefinition. We could also write a more general theory by modifying the scalar field Lagrangian (79 1 to contain 
non-canonical kinetic terms. 



(82) 



where F{) is some function of the field and its derivatives. Such Lagrangians appear frequently in models of inflation 
based on string theory, and are a topic of considerable current research interest. We could also complicate the 
gravitational sector by replacing the Ricci scalar R with a more complicated function / (R). An example of such a 
model is the inflation model of Starobinsky [5T], which can be reduced to the form (78 1 through a field redefinition. 
We could also introduce multiple scalar fields. 



Here we will confine ourselves for simplicity to a canonical Lagrangian ( 79 1 of a single scalar field, for which the 
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only adjustable quantity is the choice of potential V {(/)). For simplicity, we assume a flat spacetime, 

/l \ 



\ 



(83) 



and the equation of motion for the field with a Lagrangian given by Eq. ( 79 1 is: 

5V 



= 0, 



(84) 



where an overdot indicates a derivative with respect to the coordinate time t, and H — a/ a is the Hubble parameter. 
We will be particularly interested in the homogeneous mode of the field, for which the gradient term vanishes, Vcfi — 0, 
so that the the functional derivative dV/Scf) simplifies to an ordinary derivative, and the equation of motion simplifies 
to4 

4> + 3H,i} + V' {(j)) = 0. 

The stress-energy for a scalar field is given by 

and, for a homogeneous field, it takes the form of a perfect fluid with energy density p and pressure p, with 

1 



(85) 
(86) 



P 
P 



:r + ^(0), 



1 • 



(87) 



We see that the de Sitter limit, p ~ — p, is just the limit in which the potential energy of the field dominates the 
kinetic energy, (j)'^ <^ V {(p). This limit is referred to as slow roll, and under such conditions the universe expands 
quasi-exponentially. 



Hdt 1 = e"^, 



a (t) cc exp 

where it is conventional to define the number of e-folds N with the sign convention 

dN = -Hdt, 

so that A'^ is large in the far past and decreases as we go forward in time and as the scale factor a increases. 
This can be made quantitative by plugging the energy and pressure (87) into the Friedmann Equation 



(89) 



3r7ipj 



and the Raychaudhuri Equation, which we write in the convenient form 

47r 



{p + 3p) = H'{l-e). 



(90) 



(91) 



* The astute reader may well ask: if we are claiming inflation is a solution to the problems of flatness and homogeneity in the universe, 
why are we assuming flatness and homogeneity from the outset? The answer is that, as long as inflation gets started somehow and 
goes on for long enough, the late-time behavior of the field <f> will always be described by Eq. \S5\ . We will see later that we only 
have observational access to the end of the inflationary period, and therefore a consistent theory of initial conditions is not required for 
investigating the observational consequences of inflation. 
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Here H is given in terms of (j) by the Friedmann Equation ( 90 1 , and the parameter e specifies the equation of state, 

\ 2 



1 = 



47r 



H 



It is a straightforward exercise to show that e is related to the evolution of the Hubble parameter by 

dlnH _ I dH 



(92) 



(93) 



where N is the number of e- folds ( 89 1 . This is a useful parameterization because the condition for accelerated 

-p is equivalent to e 



expansion a > is simply equivalent to e < 1 . The de Sitter limit p 
V (0) dominates the energy density, and 



0, so that the potential 



3m|[ 



We make the additional approximation that the friction term in the equation of motion ( 85 ) dominates, 

so that the equation of motion for the scalar field is approximately 

3i7</> + V (0) ~ 0. 



(94) 



(95) 



(96) 



Equation (96) together with the Friedmann Equation (94 1 are together referred to as the slow roll approximation. 
The condition (95) can be expressed in terms of a second dimcnsionless parameter, conventionally defined as 



H(t) 



1 de 
2edN' 



(97) 



The parameters e and 77 are referred to as slow roll parameters, and the slow roll approximation is valid as long as 
both are small, e, It^I ^ 1. It is not obvious that this will be a valid approximation for situations of physical interest: 
77 need not be small for infiation to take place. Inflation takes place when e < 1, regardless of the value of rj. We later 
demonstrate explicitly that slow roll does in fact hold for interesting choices of inflationary potential. In the limit of 
slow roll, we can use Eqs. (94 96 1 to write the parameter e approximately as 



47r 
mil 



H 



'"■PI 
167r 



y(0) 



(98) 



The inflationary limit, e <C 1 is then just equivalent to a fleld evolving on a flat potential, V' {(f)) ^ V {(f)). The second 
slow roll parameter 77 can likewise be written approximately as: 



V = 



H(f> 



'"■PI 
87r 



V{^) 



V'{4>) 
V{(f>) 



(99) 



so that the curvature V" of the potential must also be small for slow roll to be a valid approximation. Similarly, we 
can write number of e-folds as a function N {(f)) of the field as: 



TV 



Hdt ■ 



H 



'PI "'0, 



V (0)' 



(100) 



The limits on the last integral are defined such that (f)^ is a fixed field value, which we will later take to be the end 
of infiation, and N increases as we go backward in time, representing the number of e-folds of expansion which take 
place between field value cf) and (f)e- 
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FIG. 14: A schematic of the potential for inflation. Inflation takes place on the region of the potential which is sufficiently 
"flat" , and reheating takes place near the true vacuum for the field. 



The qualitative picture of scalar field-driven inflation is that of a phase transition with order parameter given by 
the field (f). At early times, the energy density of the universe is dominated by the field (j) which is slowly evolving 



on a nearly constant potential, so that it approximates a cosmological constant (Fig. 14 1. During this period, the 
universe is exponentially driven toward flatness and homogeneity. Inflation ends as the potential steepens and the 
field begins to oscillate about its vacuum state at the minimum of the potential. At this point, we have an effectively 
zero-temperature scalar in a state of coherent oscillation about the minimum of the potential, and the universe is a 
huge Bose-Einstein condensate: hardly a hot Big Bang! In order to transition to a radiation-dominated hot Big Bang 
cosmology, the energy in the infiaton field must decay into Standard Model particles, a process generically termed 
reheating. This process is model-dependent, but it typically happens very rapidly. Note that the field (j) need not 
be a fundamental field like a Higgs boson (although it could in fact be fundamental). Any order parameter for a 
phase transition will do, as long as it has the quantum numbers of vacuum, and the effective potential has the correct 
properties. The infiaton cj) could well be a scalar composite of more fundamental degrees of freedom, the coordinate 
of a brane in a higher-dimensional compactification from string theory, a supersymmetric modulus, or something 
even more exotic. The simple single-field picture we discuss here is therefore an effective representation of a large 
variety of underlying fundamental theories. All of the physics important to infiation is contained in the shape of the 
potential V {(j)). (The details of the underlying theory are important for understanding the epoch of reheating, since 
the reheating process depends crucially on the specific couplings of the infiaton to the other degrees of freedom in the 
theory.) 

How long does infiation need to go on in order to solve the fiatness and horizon problems? We use a thermodynamic 
argument, which rests on a simple fact about cosmological expansion: as long as there are no decays or annihilations of 
massive particles, all other interactions conserve photon number, so that the number of photons in a comoving volume 
is constant. Since the entropy of photons is proportional to the number density, that means the entropy per comoving 
volume is also constant. Therefore, the total entropy in the Cosmic Microwave Background (or, equivalently, the total 
number of photons) is a convenient measure of spatial volume in the universe. Since the entropy per photon s is (up 
to a few constants) given by the cube of the temperature, 

s~T3, (101) 

the total photon entropy S in our current horizon volume is of order 

5w-^aMB4-(^)'-lo«^ (102) 

where we have taken the CMB temperature to be 2.7 K and the current Hubble parameter Hq to be 70 km/s/MpC. 
(The interesting unit conversion from km/s/MpC to Kelvin is left as an exercise for the reader.) 

Let us consider a highly over-simplified picture of the universe, in which no particle decays or annihilations occur 
between the end of infiation and today. In that case, the only time when the photon number (and therefore the 
entropy) in the universe changes is during the reheating process itself, when the infiation (j) decays into radiation and 
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sets the initial state for the hot Big Bang. Therefore, we must at minimum create an entropy of 10** during reheating. 
Let us say that the energy density during inflation is 

p^y(0)~A^ (103) 

where A is some energy scale. Therefore, the horizon size during inflation is then 

dn-H-'^^, (104) 

so that the initial volume of the inflationary "patch" which undergoes exponential expansion is 

V.-d^^"^. (105) 

Suppose inflation continues for N e-folds of expansion, so that the scale factor a increases by a factor of during 
inflation. The proper volume of the initial inflationary patch increases by the cube of the scale factor 

y,^e3^4-e^^^. (106) 

Inflation takes a tiny patch of the universe and blows it up exponentially large, but in such a way that the energy 
density remains approximately constant: we have created an exponential amount of energy out of nothing! During 
reheating, this huge store of energy in the coherently oscillating field </> decays into radiation and the temperature and 
entropy of the universe undergo an explosive increase. If reheating is highly efficient, then all or most of the energy 
stored in the inflaton field will be transformed into radiation, and the temperature of the universe after reheating will 
be of order the energy density of the inflaton field, 

Trh - A. (107) 

The entropy per comoving volume after reheating will then be srh ~ Trjj ~ , and the total entropy in our inflating 
patch will be 

S^H-VfTi^^e'^^. (108) 

Since this is our only source of entropy in our toy- model universe, this entropy must be at least as large as the entropy 
in our current horizon volume, Sjiii > 10**. The only adjustable parameter is the number of e-folds of inflation. 
Taking the logarithm of both sides gives a lower bound on N, 

iV>68 + ln( — ). (109) 
\mpij 

We will see later that the amplitude of primordial density fluctuations dp/ p ~ 10^^ typically constrains the inflationary 
energy scale to be of order A ^ lO^^'mpi, so that we have a lower limit on the number of e-folds of inflation of 

N > 7V„,i„ ~ 60. (110) 

Most inflation models hugely oversaturate this bound, with A'tot S> N^in. There is in fact no upper bound on 
the number of e-folds of inflation, an idea which is central to Linde's idea of "eternal" inflation [531 [Ml in 
which inflation, once initiated, never completely ends, with reheating occurring only in isolated patches of the cosmos. 
Furthermore, it is easy to see that our oversimplifled toy model of the universe gives a remarkably accurate estimate of 
A^min- In the real universe, all sorts of particle decays and annihilations happen between the end of inflation and today. 



which create additional entropy. However, our lower bound (1091 is only logarithmically sensitive to these processes. 
The dominant uncertainty is in the reheat temperature: it is possible that the energy scale of inflation is very low, or 
that the reheating process is very inefficient, and there are very few observational howads on these scales. We do know 
that the universe has to be radiation dominated and in equilibrium by the time primordial nucleosynthesis happens 
at temperatures of order MeV. Furthermore, the baryon asymmetry of the universe is at least a good hint that the 
Big Bang was hot to at least the scale of electroweak unification. A typical assumption is that the reheat temperature 
is something between 1 TeV and 10^^ GeV, which translates into a range for A^min of order [5S1 [57] 

A^„,i„ ~ [46, 60] . (Ill) 
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A. Example: the potential 

We are now in a position to apply this to a specific case. We use the simple case of a quartic potential, 

V{<j)) = Xcl)'^. (112) 



The slow roll equations (96 94 1 imply that the field evolves as: 



V (0) ^ /mg.1 
3H V 247r yvT^y 



PI ■ cx</>. (113) 



Note that this potential does not much qualitatively resemble the schematic in Fig. 13 the "flatness" of the potential 



arises because the energy density V {4>) ex (j) rises much more quickly than the kinetic energy, (jj cx (jr, so that if 
the field is far enough out on the potential, the slow roll approximation is self-consistent. The field rolls down to the 
potential toward the vacuum at the origin, and the equation of state is determined by the parameter e. 

The field value (jj^ at the end of inflation is when e {(pe) = 1, or 



(115) 



For (j) > (j)e, e < 1 and the universe is inflating, and for < 0e, e > 1 and the expansion enters a decelerating phase. 
Therefore, even this simple potential has the necessary characteristics to support a period of early-universe inflation 
followed by reheating and a hot Big Bang cosmology. What about the requirement that the universe inflate for at 



least 60 e-folds? Using Eq. (98 1, we can express the number of e-folds before the end of inflation (100 1 as 



N=^^ ^= = 7r — -1, 116 



where we integrate backward from 0e to (f> to be consistent with the sign convention (89 1. Therefore the field value N 
e-folds before the end of inflation is 

(j)N = mpi^J^^-^, (117) 

so that 

(/)60 = 4.4TOP1. (118) 

We obtain sufficient inflation, but at a price: the field must be a long way (several times the Planck scale) out on 
the potential. However, we do not necessarily have to invoke quantum gravity, since for small enough coupling A, the 
energy density in the field can much less than the Planck density, and the energy density is the physically important 
quantity. 

In this section, we have seen that the basic picture of an early epoch in the universe dominated by vacuum-like 
energy, leading to nearly exponential expansion, can be realized within the context of a simple scalar field theory. The 
equation of state for the field approximates a cosmological constant p = — p when the energy density is dominated 
by the field potential V {(j)), and inflation ends when the potential becomes steep enough that the kinetic energy 
(p^ /2 dominates over the potential. To solve the horizon and flatness problems and create a universe consistent with 
observation, we must have at least 60 or so e-folds of inflation, although in principle inflation could continue for much 
longer than this minimum amount. This dynamical explanation for the flatness and homogeneity of the universe is an 
interesting, but hardly compelling scenario. It could be that the universe started out homogeneous and flat because 
of initial conditions, either through some symmetry we do not yet understand, or because there are many universes, 
and we just happen to flnd ourselves in a highly unlikely realization which is homogeneous and geometrically flat. 
In the absence of any other observational handles on the physics of the very early universe, it is impossible to tell. 
However, flatness and homogeneity are not the whole story: inflation provides an elegant mechanism for explaining 
the inhomogeneity of the universe as well, which we discuss in Section |V) 
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V. PERTURBATIONS IN INFLATION 

The universe we live in today is homogeneous, but only when averaged over very large scales. On small scales, the 
size of people or solar systems or galaxies or even clusters of galaxies, the universe we see is highly inhomogeneous. 
Our world is full of complex structure, created by gravitational instability acting on tiny "seed" perturbations in the 
early universe. If we look as far back in time as the epoch of recombination, the universe on all scales was homogeneous 
to a high degree of precision, a few parts in 10^. Recent observational efforts such as the WMAP satellite have made 
exquisitely precise measurements of the first tiny inhomogeneities in the universe, which later collapsed to form the 



structure we see today. (We discuss the WMAP observation in more detail in Section VI ) Therefore, another mystery 
of Big Bang cosmology is: what created the primordial perturbations? This mystery is compounded by the fact that 
the perturbations we observe in the CMB exhibit correlations on scales much larger than the horizon size at the 
time of recombination, which corresponds to an angular multipole of ^ ~ 100, or about 1° as observed on the sky 
today. This is another version of the horizon problem: not only is the universe homogeneous on scales larger than the 
horizon, but whatever created the primordial perturbations must also have been capable of generating fluctuations 
on scales larger than the horizon. Inflation provides just such a mechanism [7l [8t [9t flO l [TT | fT2 j [T3|. 

Consider a perturbation in the cosmological fluid with wavelength A. Since the proper wavelength redshifts with 
expansion, Aprop oc a(t), the comoving wavelength of the perturbation is a constant, Acom — const. This is true of 
photons or density perturbations or gravitational waves or any other wave propagating in the cosmological background. 
Now consider this wavelength relative to the size of the horizon: We have seen that in general the horizon as measured 
in comoving units is proportional to the conformal time, du cx r. Therefore, for matter- or radiation-dominated 
expansion, the horizon size grows in comoving units, so that a comoving length which is larger than the horizon 
at early times is smaller than the horizon at late times: modes "fall into" the horizon. The opposite is true during 
inflation, where the conformal time is negative and evolving toward zero: the comoving horizon size is still proportional 
to T, but it now shrinks with cosmological expansion, and comoving perturbations which are initially smaller than 
the horizon are "redshifted" to scales larger than the horizon at late times (Fig. 
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If the universe is inflating at early times, and radiation- or matter-dominated at late times, perturbations in the 
density of the universe which are initially smaller than the horizon are redshifted during inflation to superhorizon scales. 
Later, as the horizon begins to grow in comoving coordinates, the perturbations fall back into the horizon, where they 
act as a source for structure formation. In this way inflation explains the observed properties of perturbations in the 
universe, which exist at both super- and sub-horizon scales at the time of recombination. Furthermore, an important 
consequence of this process is that the last perturbations to exit the horizon are the first to fall back in. Therefore, 
the shortest wavelength perturbations are the ones which exited the horizon just at the end of inflation, = 0, and 
longer wavelength perturbations exited the horizon earlier. Perturbations about the same size as our horizon today 
exited the horizon during inflation at around A^ = 60. Perturbations which exited the horizon earlier than that, 
N > 60, are still larger than our horizon today. Therefore, it is only possible to place observational constraints on 
the end of inflation, about the last 60 e-folds. Everything that happened before that, including anything that might 
tell us about the initial conditions which led to inflation, is most probably inaccessible to us. 

This kinematic picture, however, does not itself explain the physical origin of the perturbations. Inflation driven by 
a scalar fleld provides a natural explanation for this as well. The inflaton field 4> evolving on the potential V ((/)) will 
not evolve completely classically, but will also be subject to small quantum fluctuations about its classical trajectory, 
which will in general be inhomogeneous. Since the energy density of the universe during inflation is dominated by 
the inflaton field, quantum fiuctuations in couple to the spacetime curvature and result in fiuctuations in the 
density of the universe. Therefore, in the same way that the classical behavior of the field provides a description 
of the background evolution of the universe, the quantum behavior of provides a description of the inhomogeneous 
perturbations about that background. We defer a full treatment of infiaton perturbations to Appendix A, and in 
the next section focus on the much simpler case of quantizing a decoupled scalar ip in an infiationary spacetime. In 
addition to its relative simplicity, this case has direct relevance to the generation of gravitational waves in inflation. 

A. The Klein-Gordon Equation in Curved Spacetime 

Consider an arbitrary free scalar fleld, which we denote if to distinguish it from the inflaton fleld (f>. The Lagrangian 
for the field is 



C = -g^'^d^ipd.ip, (119) 
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FIG. 15: A conformal diagram of the horizon in an inflationary universe. The comoving horizon shrinks during inflation, 
and grows during the radiation- and matter-dominated expansion, while the comoving wavelengths of perturbations remain 
constant. This drives comoving perturbations to "superhorizon" scales. 



and varying the action (1191 gives the Euler-Lagrange equation of motion 

1 



V=5 



(120) 



It will prove convenient to express the background FRW metric in conformal coordinates 

g^iu = (t) ■q^.u (121) 

instead of the coordinate-time metric ( 83 1 we used in Section ( |IV[ ) . Here r is the conformal time and rj^y — 

raina 



diag. (1, — 1, — 1, — 1) is the Minkowski metric. In conformal coordinates, the free scalar equation of motion (1201 
is 



2( - )^'-VV = 0, 



(122) 



where ' = d/dr is a derivative with respect to conformal time. Note that unlike the case of the inflaton 0, we are 
solving for perturbations and therefore retain the gradient term V^tp. The field is a decoupled spectator field 
evolving in a fixed cosmological background, and does not effect the time evolution of the scale factor a(T). An 
example of such a field is gravitational waves. If we express the spacetime metric as an FRW background plus 
perturbation 5g^u, we can express the tensorial portion of the perturbation in general as a sum of two scalar degrees 
of freedom 



5gai = SgiQ = 
327r 



(123) 
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where i,j = 1,2,3, and e^'^ are longitudinal and transverse polarization tensors, respectively. It is left as an exercise 



for the reader to show that the scalars (p+^x behave to linear order as free scalars, with equation of motion (1221 



To solve the equation of motion ( 122 1, we first Fourier expand the field into momentum states (pk 

W^. (r) feke^'^-- + (r) 6^e-*-] . (124) 

Note that the coordinates x are comoving coordinates, and the wavevector k is a comoving wavevector, which does 
not redshift with expansion. The proper wavevector is 

kprop = k/a (r) . (125) 

Therefore, the comoving wavenumber k is not itself dynamical, but is just a set of constants labeling a particular 
Fourier component . The equation of motion for a single mode Lp\^ is 

(^;^ + 2 (^^^ (^;, + fcVk = o. (126) 

It is convenient to introduce a field redefinition 

life EE a (r) (^k (t) , (127) 
and the mode function Uk obeys a generalization of the Klein-Gordon equation to an expanding spacetime, 

a" 



k' 

a 



Uk - 0. (128) 



(We have dropped the vector notation k on the subscript, since the Klein-Gordon equation depends only on the 
magnitude of fc.) 

Any mode with a fixed comoving wavenumber k redshifts with time, so that early time corresponds to short 
wavelength (ultraviolet) and late time corresponds to long wavelength (infrared). The solutions to the mode equation 
show qualitatively different behaviors in the ultraviolet and infrared limits: 

• Short wavelength limit, k 3> a" /a. In this case, the equation of motion is that for a conformally Minkowski 
Klein- Gordon field. 



with solution 



< -I- k^Uk = 0, (129) 

ut (r) = ^ (Afce-^'^^ + Bke'''^) . (130) 

Note that this is in terms of conformal time and comoving wavenumber, and can only be identified with an 
exactly Minkowski spacetime in the ultraviolet limit. 

• Long wavelength limit, k <^ a" /a. In the infrared limit, the mode equation becomes 

a"uk — au'l, (131) 

with the trivial solution 

Uk oc a ^ ipk = const. (132) 

This illustrates the phenomenon of mode freezing: field modes ipk with wavelength longer than the horizon size 
cease to be dynamical, and asymptote to a constant, nonzero amplitude.^ This is a quantitative expression of 
our earlier qualitative notion of particle creation at the cosmological horizon. The amplitude of the field at long 
wavelength is determined by the boundary condition on the mode, i.e. the integration constants Ak and Bk. 

Therefore, all of the physics boils down to the question of how we set the boundary condition on field perturbations in 
the ultraviolet limit. This is fortunate, since in that limit the field theory describing the modes becomes approximately 
Minkowskian, and we know how to quantize fields in Minkowski Space. Once the integration constants are fixed, the 
behavior of the mode function Uk is completely determined, and the long-wavelength amplitude of the perturbation 
can then be calculated without ambiguity. We next discuss quantization. 



The second solution to this equation is a decaying mode, which is always subdominant in the infrared limit. 
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B. Quantization 



We have seen that the equation of motion for field perturbations approaches the usual Minkowski Space Klein- 
Gordon equation in the ultraviolet limit, which corresponds to the limit of early time for a mode redshifting with 
expansion. We determine the boundary conditions for the mode function via canonical quantization. To quantize 



the field ipk, we promote the Fourier coefficients in the classical mode expansion (124) to annihilation and creation 
operators 



l.'i 



with commutation relation 



6k, 



= 6^ (k - k') 



(133) 



(134) 



Note that the commutator in an FRW background is given in terms of comoving wavenumber, and holds whether we 
are in the short wavelength limit or not. In the short wavelength limit, this becomes equivalent to a Minkowski Space 
commutator. The quantum field ip is then given by the usual expansion in operators 6k, 



(/3(t,x) = 

The corresponding canonical momentum is 



{2n) 



3/2 



[V'k (r) b. 



ke 



H.C. 



n(T,x) = 



sc 



dtp 



6 (doip) dr ' 

It is left as an exercise for the reader to show that the canonical commutation relation 

[^(r,x),n(r,x')]-*<5-Mx-x') 
corresponds to a Wronskian condition on the mode uj-, 



Uk 



dr 



= dr 



which for the ultraviolet mode function ( 130 1 results in a condition on the integration constants 

|2 



\Bk 



1. 



(135) 



(136) 



(137) 



(138) 



(139) 



This quantization condition corresponds to one of the two boundary conditions which are necessary to completely 
determine the solution. The second boundary condition comes from vacuum selection, i.e. our definition of which 
state corresponds to a zero-particle state for the system. In the next section, we discuss the issue of vacuum selection 
in detail. 



C. Vacuum Selection 



Consider a quantum field in Minkowski Space. The state space for a quantum field theory is a set of states 
|n(ki), . . . , n(k.i)) representing the number of particles with momenta ki,...,kj. The creation and annihilation 
operators and fik act on these states by adding or subtracting a particle from the state: 

4l"(k)) -\^n|n(k) + l) 

flk |n(k)) = |n(k) - 1) . (140) 
The ground state, or vacuum state of the space, is just the zero particle state: 

ak|0)-0. (141) 
Note in particular that the vacuum state |0) is not equivalent to zero. The vacuum is not nothing: 

|0) ^ 0. (142) 
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To construct a quantum field, we look at the familiar classical wave equation for a scalar field, 

V^^^ 0. 

To solve this equation, we decompose into Fourier modes Wk, 



where the mode functions Uk(t) satisfy the ordinary differential equation 

itk + fc^Uk = 0. 



(143) 



(144) 



(145) 



This is a classical wave equation with a classical solution, and the Fourier coefficients Ok are just complex numbers. 
The solution for the mode function is 



uic oc e 



-iuJkt 



where ujk satisfies the dispersion relation 



Lol-k^^ 0. 



To turn this into a quantum field, we identify the Fourier coefficients with creation and annihilation operators 



and enforce the commutation relations 



^S^k- k') 



(146) 
(147) 
(148) 
(149) 



This is the standard quantization of a scalar field in Minkowski Space, which should be familiar. But what probably 
is not familiar is that this solution has an interesting symmetry. Suppose we define a new mode function which is 
a rotation of the solution (146): 



Uk = A{k)e 



B{k)e 



iujt — i\i.--x. 



(150) 



This is also a perfectly valid solution to the original wave equation ( 143 1, since it is just a superposition of the Fourier 



modes. But we can then re-write the quantum field in terms of our original Fourier modes and new operators fek and 



bt and the original Fourier modes e*''" as: 



where the new operators 6k are given in terms of the old operators Ok by 

fek = A(fc)ak + S*(fc)a[. 



(151) 



(152) 



This is completely equivalent to our original solution ( |144[ ) as long as the new operators satisfy the same commutation 
relation as the original operators. 



(k - k') 



This can be shown to place a condition on the coefficients A and B, 

\Af-\Bf = l. 



(153) 



(154) 



Otherwise, we are free to choose A and B as we please. 

This is just a standard property of linear differential equations: any linear combination of solutions is itself a 
solution. But what does it mean physically? In one case, we have an annihilation operator Ok which gives zero when 
acting on a particular state which we call the vacuum state: 



akiOa) = 0. 



(155) 
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Similarly, our rotated operator gives zero when acting on some state 

&k|Ob)=0. (156) 
The point is that the two "vacuum" states are not the same 

|0a> ^ \0b) . (157) 

From this point of view, we can define any state we wish to be the "vacuum" and build a completely consistent 
quantum field theory based on this assumption. From another equally valid point of view this state will contain 
particles. How do we tell which is the physical vacuum state? To define the real vacuum, we have to consider the 
spacetime the field is living in. For example, in regular special relativistic quantum field theory, the "true" vacuum 
is the zero-particle state as seen by an inertial observer. Another more formal way to state this is that we require the 
vacuum to be Lorentz symmetric. This fixes our choice of vacuum |0) and defines unambiguously our set of creation 
and annihilation operators a and . A consequence of this is that an accelerated observer in the Minkowski vacuum 
will think that the space is full of particles, a phenomenon known as the Unruh effect [5S]. The zero-particle state for 
an accelerated observer is different than for an inertial observer. The case of an FRW spacetime is exactly analogous, 
except that the FRW equivalent of an inertial observer is an observer at rest in comoving coordinates. Since an FRW 
spacetime is asymptotically Minkowski in the ultraviolet limit, we choose the vacuum field which corresponds to the 
usual Minkowski vacuum in that limit, 

Uk (r) cx 6-''=" ^Ak = l, Bk = 0. (158) 

This is known as the Bunch-Davies vacuum. This is not the only possible choice, although it is widely believed to 
be the most natural. The issue of vacuum ambiguity of infiationary perturbations is a subject which is extensively 
discussed in the literature, and is still the subject of controversy. It is known that the choice of vacuum is potentially 
sensitive to quantum-gravitational physics jSHl ISni EI], a subject which is referred to as Trans-Planckian physics 
[m [621 [63] . For the remainder of our discussion, we will assume a Bunch-Davies vacuum. 

The key point is that quantization and vacuum selection together completely specify the mode function, up to an 
overall phase. This means that the amplitude of the mode once it has redshifted to long wavelength and frozen 
out is similarly determined. In the next section, we solve the mode equation at long wavelength for an infiationary 
background. 



D. Exact Solutions and the Primordial Power Spectrum 



The exact form of the solution to Eq. (1281 depends on the evolution of the background spacetime, as encoded in 



a (t), which in turn depends on the equation of state of the field driving inflation. We will consider the case where the 
equation of state is constant, which will not be the case in general for scalar field-driven inflation, but will nonetheless 
turn out to be a good approximation in the limit of a slowly rolling fleld. Generalizing Eq. (77 1 to the case of arbitrary 
equation of state parameter e = const., the conformal time can be written 



1 



1 



1 - e 



and the Fricdmann and Raychaudhuri Equations ( 14 1 give 

a^H^ (2 - e) , 



a 

a 



(159) 



(160) 



where a prime denotes a derivative with respect to conformal time. The conformal time, as in the case of de Sitter 
space, is negative and tending toward zero during inflation. (Proof of these relations is left as an exercise for the 
reader.) We can then write the mode equation (1281 as 



[fc2 - a^H^ (2 - e)] Uk = 0. 



(161) 



Using Eq. (159) to write aH in terms of the conformal time r, the equation of motion becomes 



T2(l-e)'<+ (fcT)'(l-e)'-(2-e) Uk = 



(162) 
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This is a Bessel equation, with solution 



where the index v is given by: 



3-e 
2(1^- 



The quantity — fcr has special physical significance, since from Eq. (1591 we can write 

k 



i-kr) (1 - e) = 



(163) 



(164) 



(165) 



where the quantity (k/aH) expresses the wavenumber k in units of the conioving horizon size dn ^ (aH)^^. Therefore, 
the short wavelength limit is —kr—^ ~oo, or (k/aH) 3> 1. The long- wavelength limit is — fcr — > 0, or (k/aH) <C 1. 

The simple case of de Sitter space (p = — p) corresponds to the limit e = 0, so that the Bessel index is v = 3/2 and 
the mode function ( 163 1 simplifies to 



Uk oc 



ki 



kr 



In the short wavelength limit, (—kr) — > —oo, the mode function is given, as expected, by 

Selecting the Bunch-Davies vacuum gives Uk oc e*'^'^, and canonical quantization fixes the normalization. 



1 



Uk 



-ikr 



2k 



Therefore, the fully normalized exact solution is 



Uk = 



1 fk 



2k V kr 



(166) 



(167) 



(168) 



(169) 



This solution has no free parameters aside from an overall phase, and is valid at all wavelengths, including after the 
mode has been redshifted outside of the horizon and becomes non-dynamical, or "frozen" . In the long wavelength 
limit, —kr 0, the mode function (169 1 becomes 



1 



Uk 



2k \{-kT) 



i ( aH 



2k 



oc a, 



consistent with the qualitative result (1321. Therefore the field amplitude ipk is given by 

H 



Wk\ 



V2fc3/^ 



= const. 



(170) 



(171) 



The quantum mode therefore displays the freezeout behavior we noted qualitatively above (Fig. 16). The amplitude 
of quantum fiuctuations is conventionally expressed in terms of the two-point correlation function of the field ip. It is 
left as an exercise for the reader to show that the vacuum two-point correlation function is given by 



(0|^(r,x)^(T,x')|0) = 



d^k 



{2n) 
dk 



where the power spectrum P (k) is defined as 

F(fc)^ 



k^ 
2^ 



Uk_ 

a 



, -kr 0. 



(172) 



(173) 
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FIG. 16: The normahzed mode function in de Sitter space, showing oscillatory behavior on subhorizon scales k/aH > 1, and 
mode freezing on superhorizon scales, k/aH < 1. 



The power per logarithmic interval k in the field fluctuation is then given in the long wavelength limit by the Hubble 
parameter H = const. This property of scale invariance is exact in the de Sitter limit. 

In a more general model, the spacetime is only approximately de Sitter, and we expect that the power spectrum of 
field fluctuations will only be approximately scale invariant. It is convenient to express this dynamics in terms of the 
equation of state parameter e, 

1 ,,,,, 

We must have e < 1 for inflation, and for a slowly rolling field jT^j ^ 1 means that e will also be slowly varying, 
e ~ const. It is straightforward to show that for e = const. ^ that: 

• The Bunch-Davies vacuum corresponds to the positive mode of Eq. ( |163 1, 

Uk oc V-fcr [J^ (-kr) + lY^ (-fcr)] . (175) 

• Quantization fixes the normalization as 



• The power spectrum in the long- wavelength limit k/aH — > is a power law in k: 

r(3/2) ^ V2vry \aH{l 



[P(k)f'^r-y^^il-e)(^)(^^Y'~\ (177) 
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where F (z/) is a gamma function, and 



2{l-ey 



(178) 



Proof is left as an exercise for the reader.^ Note that in the case e — const., both the background and perturbation 
equations are exactly solvable. 

We can use these solutions as approximate solutions in the more general slow roll case, where e <C 1 ~ const., so 
that the dependence of the power spectrum on k is approximately a power-law. 



P{k) cx fc". 



(179) 



with spectral index 



n — i — 2v = i — 



1 -e 



-2e. 



(180) 



Equation (1771 is curious, however, because it does not obviously exhibit complete mode freezing at long wavelength, 
since a and H both depend on time. We can show that P (k) does in fact approach a time-dependent value at long 
wavelength by evaluating 



_d_ 
dN 



H 



k 



3/2-1/' 



H 



= He 

-0, 



aH 



_d_ 
dN 



1-e \ aH 



aH 



k ek 



(181) 



which can be easily shown using a oc and H (x e*^ . That is, the time-dependent quantities a and H in Eq. (1771 
are combined in such a way as to form an exactly conserved quantity. Since it is conserved, we are free to evaluate it 
at any time (or value of aH) that we wish. It is conventional to evaluate the power spectrum at horizon crossing, or 
at aH = fc, so that 



P"' (k) 



27T J 



(182) 



k=aH 



where we have approximated the z/-dependent multiplicative factor as order unity. ^ 

It is straightforward to calculate the spectral index ( 180 1 directly from the horizon crossing expression ( 182 1 by 
using 



a (X e ^ , H (X e'^ 

so that we can write derivatives in k at horizon crossing as derivatives in the number of e-folds TV, 

1 d{aH) 



dlnk\j,=aH = dln{aH) = 
The spectral index is then, to lowest order in slow roll 

dlnP{k) 

n = 

din fc 



aH dN 



-dN = (e - 1) dN. 



(183) 



(184) 



fc dH^ 

milk 
1 



k=aH 

dH^ 



H^{e-l) dN 



® Note that the quantization condition ( |l37| can be apphed to the solution | |l63| exactly, resulting in the normalization condition l |l39| , 
without approximating the solution in the short- wavelength limit! 

^ This is not the value of the scalar field power spectrum at the moment the mode is physically crossing outside the horizon, as is often 
stated in the literature: it is the value of the power spectrum in the asymptotic long-wavelength limit. It is easy to show from the exact 
solution | |176| l that the mode function is still evolving with time as it crosses the horizon at A: = aH, and the asymptotic amplitude 
differs from the amplitude at horizon crossing by about a factor of two. See Ref. 1641 for a more detailed discussion of this point. 
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2e 

-2e, 



(185) 



in agreement with ( |180[ ). Note that we are rather freely changing variables from the wavenumber k to the comoving 
horizon size (aH)^ to the number of e-folds N. As long as the cosmological evolution is monotonic, these are all 
different ways of measuring time: the time when a mode with wavenumber k exits the horizon, the time at which the 
horizon is a particular size, the number of e-folds N and the field value are all effectively just different choices of 
a clock, and we can switch from one to another as is convenient. For example, in the slow roll approximation, the 
Hubble parameter H is just a function of (j), H cc y^V (0). Because of this, it is convenient to define N (fc) to be the 
number of e-folds ( 100 1 when a mode with wavenumber k crosses outside the horizon, and 0jv (k) to be the field value 
N (fc) e-folds before the end of inflation. Then the power spectrum can be written equivalently as either a function of 
k or of (p: 
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(186) 



Wavenumbers k are conventionally normalized in units of /iMpc~ as measured in the current universe. We can 
relate N to scales in the current universe by recalling that modes which are of order the horizon size in the universe 
today, k ~ ciqHq, exited the horizon during inflation when N = [46,60], so that we can calculate the amplitude of 
perturbations at the scale of the CMB quadrupole today by evaluating the power spectrum for field values between 
046 and 060. 

One example of a free scalar in inflation is gravitational wave modes, where the transverse and longitudinal po- 
larization states of the gravity waves evolve as independent scalar fields. Using Eq. (123), we can then calculate 



the power spectrum in gravity waves (or tensors) as the sum of the two-point correlation functions for the separate 
polarizations: 



Pt = (<54>=2 
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with spectral index 



ut = — 2e. 



(187) 
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If the amplitude is large enough, such a spectrum of primordial gravity waves will be observable in the cosmic 
microwave background anisotropy and polarization, or be directly detectable by proposed experiments such as Big 
Bang Observer [65l |66] . 

The second type of perturbation generated during inflation is perturbations in the density of the universe, which are 
the dominant component of the CMB anisotropy ST/T Sp/p ~ 10~^, and are responsible for structure formation. 
Density, or scalar perturbations are more complicated than tensor perturbations because they are generated by 
quantum fluctuations in the inflaton field itself: since the background energy density is dominated by the inflaton, 
fluctuations of the inflaton up or down the potential generate perturbations in the density. The full calculation requires 
self-consistent General Relativistic perturbation theory, and is presented in Appendix A. Here we simply state the 
result: Perturbations in the inflaton field S4> — H/2t: generate density perturbations with power spectrum 
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k=aH 



where N is the number of e-folds. Scalar perturbations are therefore enhanced relative to tensor perturbations by a 
factor of 1/e. The scalar power spectrum is also an approximate power-law, with spectral index 



ns - I = 



{e-\)dN\ e 
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(190) 



where rj is the second slow roll parameter (97 1. Therefore, for any particular choice of inflationary potential, we have 



four measurable quantities: the amplitudes Pt and P-ji of the tensor and scalar power spectra, and their spectral 
indices and ng. However, not all of these parameters are independent. In particular, the ratio r between the scalar 
and tensor amplitudes is given by the parameter e, as is the tensor spectral index ut'. 
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This relation is known as the consistency condition for single-field slow roll inflation, and is in principle testable by a 
sufliciently accurate measurement of the primordial perturbation spectra. 

In the next section, we apply these results to our example A^'' potential and calculate the inflationary power spectra. 



E. Example: X(fi'^ 



For the case of our example model with V {(f)) — A^"*, it is now straightforward to calculate the scalar and tensor 
perturbation spectra. We express the normalization of the power spectra as a function of the number of e-folds N by 
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For perturbations about 



where we have used the slow roll expressions for _ff (94l and e (98l and Eq. (117) for 
the current size of our horizon, N ~ 60, and CMB normalization forces the self-coupling to be very small, 

A - 10"^^ 



(193) 



The presence of an extremely small parameter is not peculiar to the X(/)* model, but is generic, and is referred to as 
the fine tuning problem for inflation. 

We can similarly calculate the tensor amplitude 
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which is usually expressed in terms of the tensor/scalar ratio 
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where we have again taken N = 60. For this particular model, the power in gravitational waves is large, about 
a quarter of the power in scalar perturbations. This is not generic, but is quite model-dependent. Some choices of 
potential predict large tensor contributions (where "large" means of order 10%), and other choices of potential predict 
very tiny tensor contributions, well below 1%. 

The tensor spectral index ut is fixed by the consistency condition (1911, but the scalar spectral index ns is an 
independent parameter because of its dependence on rj: 
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The spectral index is then 



n = 1 - — ~ 0.95. (199) 



Note that we have assumed slow roll from the beginning in the calculation without a priori knowing that it is a good 
approximation for this choice of potential. However, at the end of the day it is clear that the slow roll ansatz was a 
good one, since e and 77 are both of order 0.01. 

Finally, we note that the energy density during inflation is characterized by a mass scale 

pi/4 _ A - Ai/^mpi - 10^5 GeV, (200) 

about the scale for which we expect Grand Unification to be important. This interesting coincidence suggests that 
the physics of inflation may be found in Grand Unified Theories (GUTs). Different choices of potential V {4>) will give 
different values for the amplitudes and shapes of the primordial power spectra. Since the normalization is fixed by the 
CMB to be Ptz ^ 10~^, the most useful observables for distinguishing among different potentials are the scalar/tensor 
ratio r and the scalar spectral index ng. In single-field infiationary models, the tensor spectral index is fixed by the 



consistency condition (1911, and is therefore not an independent parameter. The consistency condition can therefore 
be taken to be a prediction of single-field infiation, which is in principle verifiable by observation. In practice, this is 
very difficult, since it involves measuring not just the amplitude of the gravitational wave power spectrum, but also 
its shape. We will see in Section VI that current data place only a rough upper bound on the tensor/scalar ratio r. 



and it is highly unlikely that any near-future measurement of primordial gravitational waves will be accurate enough 
to constrain nx well enough to test the consistency condition. In the next section, we discuss current observational 
constraints on the form of the inffationary potential. 



VI. OBSERVATIONAL CONSTRAINTS 



Our simple picture of inffation generated by single, minimally coupled scalar field makes a set of very definite 
predictions for the form of primordial cosmological fiuctuations: 

• Gaussianity: since the two-point correlation function of a free scalar field (<p^) is Gaussian, cosmological pertur- 
bations generated in a single-field inflation model will by necessity also form a Gaussian random distribution. 

• Adiabaticity: since there is only one order parameter (j) governing the generation of density perturbations, we 
expect the perturbations in all the components of the cosmological fluid (baryons, dark matter, neutrinos) to 
be in phase with each other. Such a case is called adiabatic. If one or more components fluctuates out of phase 
with others, these are referred to as isocurvature modes. Single-field inflation predicts an absence of isocurvature 
fluctuations. 

• Scale invariance: In the limit of de Sitter space, fluctuations in any quantum fleld are exactly scale invariant, 
n = 1, as a result of the fact that the Hubble parameter is exactly constant. Since slow-roll inflation is quasi- 
de Sitter, we expect the perturbation spectra to be nearly, but not exactly scale invariant, with jri^ — 1| = 
|2r7-4e| < 1. 

• Scalar perturbations dominate over tensor perturbations, r — 16e. 

Furthermore, given a potential V (0), we have a "recipe" for calculating the form of the primordial power spectra 
generated during inflation: 

1. Calculate the fleld value at the end of inflation 0e from 



16t: \V{(j)e) 

2. Calculate the field value N e-folds before the end of inflation (p^ by integrating backward on the potential from 

4' = 4'e, 

TV = / , ^ . (202 
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3. Calculate the normalization of the scalar power spectrum by 



1/2 



H 



- 10"^ (203) 

<P=4>N 



where the CMB quadrupole corresponds to roughly A'' = 60. A more accurate calculation includes the uncer- 
tainty in the reheat temperature, which gives a range N ~ [46,60], and a corresponding uncertainty in the 
observable parameters. 

4. Calculate the tensor/scalar ratio r and scalar spectral index ns at N — [46, 60] by 

r= 16e((/.jv), (204) 

and 

= l-4e(0A,) + 2?7(0iv), (205) 
where the second slow roll parameter 77 is given by: 



an 



(206) 



The key point is that the scalar power spectrum P-ji and the tensor power spectrum Pt are both completely determined 
by the choice of potential V {(/)). ^ Therefore, if we measure the primordial perturbations in the universe accurately 
enough, we can in principle constrain the form of the inflationary potential. This is extremely exciting, because it 
gives us a very rare window into physics at extremely high energy, perhaps as high as the GUT scale or higher, far 
beyond the reach of accelerator experiments such as the Large Hadron Collider. 

It is convenient to divide the set of possible single-field potentials into a few basic types [55] : 



Large-field potentials (Fig. 17 1. These are the simplest potentials one might imagine, with potentials of the 
form V {(/)) — m'^(p'^, or our example case, V {(j)) — Xcj)*. Another widely-noted example of this type of model 
is inflation on an exponential potential, V {(j>) — A^ exp ((/)//i), which has the useful property that both the 
background evolution and the perturbation equations are exactly solvable. In the large-fleld case, the fleld is 
displaced from the vacuum at the origin by an amount of order ^ mp\ and rolls down the potential toward 
the origin. Large-fleld models are typically characterized by a "red" spectral index ns < I, and a substantial 
gravitational wave contribution, r ^ 0.1. 



Small-field potentials (Fig. 18 1. These are potentials characteristic of spontaneous symmetry breaking phase 
transitions, where the fleld rolls off an unstable equilibrium with V' (0) = toward a displaced vacuum. Ex- 
amples of small-fleld inflation include a simple quadratic potential, V {(p) — A (0^ — /i^) , inflation from a 
pseudo-Nambu-Goldstone boson or a shift symmetry in string theory (called Natural Inflation) with a potential 
typically of the form V {(p) = A"* [1 -I- cos {(f)/ fi)], or Coleman- Weinberg potentials, V {(p) = Xcp* In {(/)). Small-fleld 
models are characterized by a red spectral index n < 1, and a small tensor/scalar ratio, r < 0.01. 



Hybrid potentials (Fig. 19 1. A third class of models are potentials for which there is a residual vacuum energy 

when the fleld is at the minimum of the potential, for example a potential like V {(f)) = X {4>^ -\- /i^)^. In this 
case, inflation will continue forever, so additional physics is required to end inflation and initiate reheating. The 
hybrid mechanism, introduced by Linde [7D] , solves this problem by adding a second fleld coupled to the inflaton 
which is stable for (p large, but becomes unstable at a critical fleld value pc near the minimum of V {(p) . During 
inflation, however, only (p is dynamical, and these models are effectively single-fleld. Typical models of this type 
predict negligible tensor modes, r ^ 0.01 and a "blue" spectrum, ns > 1, which is disfavored by the data, and 
we will not discuss them in more detail here. (Ref. |48j contains a good discussion of current limits on general 
hybrid models.) Note also that such potentials will also support large-fleld inflation if the fleld is displaced far 
enough from its minimum. 



Strictly speaking, this is true only for scalar fields with a canonical kinetic term, where the speed of sound of perturbations is equal to 
the speed of light. More complicated scenarios such as DBI inflation |67) require specification of an extra free function, the speed of 
sound cg {</>), to calculate the power spectra. For constraints on this more general class of models, see Ref. 1681 . 
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FIG. 17: A schematic of a large-field potential. 




FIG. 18: A schematic of a small-field potential. 

An important feature of all of these models is that each is characterized by two basic parameters, the "height" 
of the potential A**, which governs the energy density during inflation, and the "width" of the potential ^. (Hybrid 
models have a third free parameter (j)c which sets the end of inflation.) In order to have a flat potential and a slowly 
rolling field, there must be a hierarchy of scales such that the width is larger than the height, A <g; /x. As we saw in 
the case of the Xcj)'^ large-field model, typical infiationary potentials have widths of order the Planck scale /i ~ mp\ 
and heights of order the scale of Grand Unification A ~ Mgut ~ 10^^ GeV, although models can be constructed for 
which infiation happens at a much lower scale [701 [ZD [Z2] • 

The quantities we are interested in for constraining models of inflation are the primordial power spectra Pji and Pt, 
which are the underlying source of the CMB temperature anisotropy and polarization. However, the observed CMB 
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FIG. 19: A schematic of a hybrid potential. 



anisotropies depend on a handful of unrelated cosmological parameters, since the primordial fluctuations are processed 
through the complicated physics of acoustic oscillations. This creates uncertainties due to parameter degeneracies: 
our best-fit values for r and 715 will depend on what values we choose for the other cosmological parameters such as 
the baryon density fib and the redshift of reionization z^i- To accurately estimate the errors on r and ns, we must fit 
all the relevant parameters simultaneously, a process which is computationally intensive, and is typically approached 
using Bayesian Monte Carlo Markov Chain techniques [73^. Here we simply show the results: Figure 20 shows the 
regions of the r, ns parameter space allowed by the WMAP 5- year data set [71] [75]. We have fit over the parameters 
f^CDM, ^^b, ^^Lambda, ^^0 1 pRi Zri, v, and Us, with a Constraint that the universe must be flat, as predicted by inflation, 
ilb + f^cDM + f^Lambda = 1- We See that the data favor a red spectrum, ns < 1, although the scale-invariant limit 
ri5 1 is still within the 95%-confidence region. Our example inflation model V (0) = Xcj)'^ is conv incingly ruled out 
by WMAP, but the simple potential V {(j)) = rn^cf)^ is nicely consistent with the data.^ Figure 
constraint with r on a logarithmic scale, with the prediction of several small-field models for reference 



shows the WMAP 
There is no 



evidence in the WMAP data for a nonzero tensor/scalar ratio r, with a 95%-confidence upper limit of r < 0.5. It is 
possible to improve these constraints somewhat by adding other data sets, for example the ACBAR high-resolution 
CMB anisotropy measurement [THI or the Sloan Digital Sky Survey [77][7H|, which improve the upper limit on the 
tensor/scalar ratio to r < 0.3 or so. Current data are completely consistent with Gaussianity and adiabaticity, as 
expected from simple single-field inflation models. In the next section, we discuss the outlook for future observation. 



VII. OUTLOOK AND CONCLUSION 



The basic hot Big Bang scenario, in which the universe arises out of a hot, dense, smooth initial state and cools 
through expansion, is now supported by a compelling set of observations, including the existence of the Cosmic 
Microwave Background, the primordial abundances of the elements, and the evolution of structure in the universe, 
all of which are being measured with unprecedented precision. However, this scenario leaves questions unanswered: 



^ Liddle and Leach point out that At/)* models are special because of their reheating properties, and should be more accurately evaluated 
at A'^ = 64 |56) . However, this assumes that the potential has no other terms which might become dominant during reheating, and in 
any case is also ruled out by WMAP5. 



WMAP 5 Limits on Inflation 
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FIG. 20: Constraints on the r, n plane from Cosmic Microwave Background measurements. Shaded regions are the regions 
allowed by the WMAP5 measurement to 68% and 95% confidence. Models plotted are "large-field" potentials V {4>) oc tj)^ and 
V (0) oc 0^ 



WMAP 5 Limits on Inflation 
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FIG. 21: Constraints on the r, n plane from Cosmic Microwave Background measurements, with the tensor/scalar ratio 
plotted on a log scale. In addition to the large-field models shown in Fig. |20[ three small-field models are plotted against the 
data: "Natural Inflation" from a pseudo-Nambu-Goldstone boson f79], with potential V {(j>) = A** [1 — cos (</)//i)], a logarithmic 
potential V {4>) oc ln{(f>) typical of supersymmetric models ^80, i81j i82j , and a Coleman- Weinberg potential V {4>) oc </)*ln((^). 
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Why is the universe so big and so old? Why is the universe so close to geometrically flat? What created the initial 
perturbations which later collapsed to form structure in the universe? The last of these questions is particularly 
interesting, because recent observations of the CMB, in particular the all-sky anisotropy map made by the landmark 
WMAP satellite, have directly measured the form of these primordial perturbations. A striking property of these 
observed primordial perturbations is that they are correlated on scales larger than the cosmological horizon at the 
time of last scattering. Such apparently acausal correlations can only be produced in a few ways |83j : 

• Inflation. 

• Extra dimensions [SI] . 

• A universe much older than Hq ^ ^55J |5B] . 

• A varying speed of light [57] . 

In addition, the WMAP data contain spectacular confirmation of the basic predictions of the inflationary paradigm: 
a geometrically flat universe with Gaussian, adiabatic, nearly scale-invariant perturbations. No other model explains 
these properties of the universe with such simplicity and economy, and much attention has been devoted to the 
implications of WMAP for inflation [li|li|Hl|73IHllHl|Mlinill3[n3[MllMllMl[lZl[13 Inflation also makes 
predictions which have not been well tested by current data but can be by future experiments, most notably a deviation 
from a scale-invariant spectrum and the production of primordial gravitational waves. A non-scale-invariant spectrum 
is weakly favored by the existing data, but constraints on primordial gravity waves are still quite poor. The outlook 
for improved data is promising: over the next flve to ten years, there will be a continuous stream of increasingly 
high-precision data made available which will allow constraint of cosmological parameters relevant for understanding 
the early universe. The most useful measurements for direct constraint of the inflationary parameter space are 
observations of the CMB, and current activity in this area is intense. The Planck satellite mission is scheduled to 
launch in 2009 [9SJ 1100] , and will be complemented by ground- and balloon-based measurements using a variety of 
technologies and strategies [43] SU [Ml [Ml |M1 IMl [M] • 

At the same time, cosmological parameter estimation is a well-developed field. A set of standard cosmological 
parameters such as the baryon density fib^^, the matter density fim^^, the expansion rate Hq = lOO/i km/sec are being 
measured with increasing accuracy. The observable quantities most meaningful for constraining models of infiation 
are the ratio r of tensor to scalar fluctuation amplitudes, and the spectral index ns of the scalar power spectrum. 
This kind of simple parameterization is at the moment sufficient to describe the highest-precision cosmological data 
sets. Furthermore, the simplest slow-roll models of inflation predict a nearly exact power-law perturbation spectrum. 
In this sense, a simple concordance cosmology is well-supported by both data and by theoretical expectation. It 
could be that the underlying universe really is that simple. However, the simplicity of concordance cosmology is at 
present as much a statement about the data as about the universe itself. Only a handful of parameters are required 
to explain existing cosmological data. Adding more parameters to the fit does no good: any small improvement 
in the fit of the model to the data is offset by the statistical penalty one pays for introducing extra parameters 
[Ml[Ml[Ml[nS[IIIl[IIl[II31[ni[IIS]- But the optimal parameter set is a moving target: as the data get better, 
we will be able to probe more parameters. It may be that a "vanilla" universe |116j of a half-dozen or so parameters 
will continue to be sufficient to explain observation. But it is reasonable to expect that, as measurements improve in 
accuracy, we will see evidence of deviation from such a lowest-order expectation. This is where the interplay between 
theory and experiment gains the most leverage, because we must understand: (1) what deviations from a simple 
universe are predicted by models, and (2) how to look for those deviations in the data. It is of course impossible to 
predict which of the many possible signals (if any) will be realized in the universe in which we live. I discuss below 
four of the best motivated possibilities, in order of the quality of current constraints. (For a more detailed treatment 
of these issues, the reader is referred to the very comprehensive CMBPol Mission Concept Study [117 .) 

Features in the density power spectrum 

Current data are consistent with a purely power-law spectrum of density perturbations, P{k) cx fc"^~^ with a "red" 
spectrum {ns < 1) favored by the data at about a 90% confidence level, a figure which depends on the choice of 
parameter set and priors. Assuming it is supported by future data, the detection of a deviation from a scale-invariant 
(ris = 1) spectrum is a significant milestone, and represents a confirmation of one of the basic predictions of infiation. 
In slow-roll inflation, this power-law scale dependence is nearly exact, and any additional scale dependence is strongly 
suppressed. Therefore, detection of a nonzero "running" a — dns/dlnk of the spectral index would be an indication 
that slow roll is a poor approximation. There is currently no evidence for scale-dependence in the spectral index, 
but constraints on the overall shape of the power spectrum are likely to improve dramatically through measurements 
of the CMB anisotropy at small angular scales, improved polarization measurements, and better mapping of large- 
scale structure. Planck is expected to measure the shape of the spectrum with 2(7 imcertainties of order An ~ 0.01 
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and Aa ~ 0.01 |118l 11191 11201 1121j . Over the longer term, measurements of 21cm radiation from neutral hydrogen 
promises to be a precise probe of the primordial power spectrum, and would improve these constraints significantly 

m- 

Primordial Gravitational Waves 

In addition to a spectrum P-jz of scalar perturbations, inflation generically predicts a spectrum Pt of tensor per- 
turbations. The relative amplitude of the two is determined by the equation of state of the fluid driving inflation, 

r = 16e (207) 

Since the scalar amplitude is known from the COBE normalization to be P-ji ^ _ff^/e ~ 10""'^°, it follows that 
measuring the tensor/scalar ratio r determines the inflationary expansion rate H and the associated energy density 
p. Typical inflation models take place with an energy density of around p ~ (10^^ GeV)*, which corresponds to a 
tensor/scalar ratio of r ~ 0.1, although this flgure is highly model-dependent. Single-fleld inflation does not make a 
definite prediction for the value of r: while many choices of potential generate a substantial tensor component, other 
choices of potential result in an unobservably small tensor/scalar ratio, and there is no particular reason to favor one 
scenario over another. 

There is at present no observational evidence for primordial gravitational waves: the current upper limit on the 
tensor/scalar ratio is around r < 0.3. Detection of even a large primordial tensor signal requires extreme sensitivity. 
The crucial observation is detection of the odd-parity, or B-mode, component of the CMB polarization signal, which 
is suppressed relative to the temperature fluctuations, themselves at the 10"'' level, by at least another four orders 
of magnitude. This signal is considerably below known foreground levels [ 123j . severely complicating data analysis. 
Despite the formidable challenges, the observational community has undertaken a broad-based effort to search for the 
B-mode, and a detection would be a boon for inflationary cosmology. Planck will be sensitive to a tensor/scalar ratio 
of around r ~ 0.1, and dedicated ground-based measurements can potentially reach limits of order r ~ 0.01. The 
proposed CMBPol polarization satellite would reach r of order 10"'^ [llTl 1124] . and direct detection experiments such 
as BBO could in principle detect r of order 10"^ [S5] . 

Primordial Non-Gaussianity 

In addition to a power-law power spectrum, inflation predicts that the primordial perturbations will be distributed 
according to Gaussian statistics. Like running of the power spectrum, non-Gaussianity is suppressed in slow-roll infla- 
tion [1251 . However, detection of even moderate non-Gaussianity is considerably more difficult. If the perturbations 
are Gaussian, the two-point correlation function completely describes the perturbations. This is not the case for non- 
Gaussian fluctuations: higher-order correlations contain additional information. However, higher-order correlations 
require more statistics and are therefore more difficult to measure, especially at large angular scales where cosmic 
variance errors are significant. Current limits are extremely weak |88l I126j . and future high angular resolution CMB 
maps will still fall well short of being sensitive to a signal from slow-roll inflation or even weakly non-slow-roll models 
|127j . It will take a strong deviation from the slow-roll scenario to generate observable non-Gaussianity. However, 
a measurement of non-Gaussianity would in one stroke rule out virtually all slow-roll inflation models and force 
consideration of more exotic scenarios such as DBI inflation [67^, Warm Inflation |128] . or curvaton scenarios |129j . 

Isocurvature perturbations 

In a universe where the matter consists of multiple components, there are two general classes of perturbation about 
a homogeneous background: adiabatic, in which the perturbations in all of the fluid components are in phase, and 
isocurvature, in which the perturbations have independent phases. Single-field infiation predicts purely adiabatic 
primordial perturbations, for the simple reason that if there is a single field responsible for infiation, then there 
is a single order parameter governing the generation of density perturbations. This is a nontrivial prediction, and 
the fact that current data are consistent with adiabatic perturbations is support for the idea of quantum generation 
of perturbations in infiation. However, current limits on the isocurvature fraction are quite weak |130( 1131] . If 
isocurvature modes are detected, it would rule out all single-field models of infiation. Multi-field models, on the 
other hand, naturally contain multiple order parameters and can generate isocurvature modes. Multi-field models 
are naturally motivated by the string "landscape" , which is believed to contain an enormous number of degrees of 
freedom. Another possible mechanism for the generation of isocurvature modes is the curvaton mechanism, in which 
cosmological perturbations are generated by a field other than the infiaton [132 , 133J. 

The rich interplay between theory and observation that characterizes cosmology today is likely to continue for the 
foreseeable future. As measurements improve, theory will need to become more precise and complete than the simple 
picture of infiation that we have outlined in these lectures, and single-field infiation models could yet prove to be a 
poor fit to the data. However, at the moment, such models provide an elegant, compelling, and (most importantly) 
scientifically useful picture of the very early universe. 
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FIG. 22: Foliations of an FRW spacetime. Comoving hypersurfaces (dashed lines) have constant density, but another choice of 
gauge (solid lines) will have unphysical density fluctuations which are an artifact of the choice of gauge. 
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APPENDIX A: THE CURVATURE PERTURBATION IN SINGLE-FIELD INFLATION 



In this section, we discuss the generation of perturbations in the density (5 (x) = 5p/p generated during inflation. 
The process is similar to the case of a free scalar field discussed in Sec. |V] the infiaton field 0, like any other scalar, 
will have quantum fluctuations which are stretched to superhorizon scales and subsequently freeze out as classical 
perturbations. The difference is that the energy density of the universe is dominated by the inflaton potential, so 
that quantum fluctuations in cj) generate perturbations in the density p. Dealing with such density perturbations 
is complicated by the fact that in General Relativity, we are free to choose any coordinate system, or gauge, we 
wish. To see why, consider the case of an FRW spacetime evolving with scale factor a (t) and uniform energy density 
p(t,x) — p{t). What we mean here by "uniform" energy density, or homogeneity, is that the density is a constant 
in comoving coordinates. But the physics is independent of coordinate system, so we could equally well work in 



coordinates t' , x' for which constant-time hypersurfaces do not have constant density (Fig. 22 1. Such a division of 
spacetime into a time coordinate and a set of orthogonal spacelike hypersurfaces is called a foliation of the spacetime, 
and is an arbitrary choice equivalent to a choice of coordinate system. 

For an FRW spacetime, comoving coordinates correspond to a foliation of the spacetime into spatial hypersurfaces 
with constant density: this is the most physically intuitive description of the spacetime. Any other choice of foliation 
of the spacetime would result in density "perturbations" which are entirely due to the choice of coordinate system. 
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Such unphysical perturbations are referred to as gauge modes. Another way to think of this is that the division 
between what we call "background" and what we call "perturbation" is itself gauge-dependent. For perturbations 
with wavelength smaller than the horizon, it is possible to define background and perturbation without ambiguity, 
since all observers can agree on a definition of time coordinate t and on an average density p (t) . Not so for super horizon 
modes: if we consider a perturbation mode with wavelength much larger than the horizon size, observers in different 
horizons will sec themselves in independently evolving, homogeneous patches of the universe: a "perturbation" can 
be defined only by comparing causally disconnected observers, and there is an inherent gauge ambiguity in how we 
do this. The canonical paper on gauge issues in General Relativistic perturbation theory is by Bardeen |134| . A good 
pedagogical treatment with a focus on inflationary perturbations can be found in Ref. |135j . 

In practice, instead of the density perturbation 6, the quantity most directly relevant to CMB physics is the 
Newtonian potential $ on the surface of last scattering. For example, this is the quantity that directly appears in 



Eq. (50 1 for the Sachs- Wolfe Effect. The Newtonian potential is related to the density perturbation S through the 



Poisson Equation: 

= AnGpa^S, (Al) 

where the factor of comes from defining the gradient V relative to comoving coordinates. Like S, the Newtonian 
potential $ is a gauge-dependent quantity: its value depends on how we foliate the spacetime. For example, we 
are free to choose spatial hypersurfaces such that the density is constant, and the Newtonian potential vanishes 
everywhere: $ {t, x) = 0. This foliation of the spacetime is equivalent to the qualitative picture above of different 
horizon volumes as independently evolving homogeneous universes. Observers in different horizons use the density p 
to synchronize their clocks with one another. Such a foliation is not very useful for computing the Sachs- Wolfe effect, 
however! Instead, we need to define a gauge which corresponds to the Newtonian limit in the present universe. To 
accomplish this, we describe the evolution of a scalar field dominated cosmology using the useful fluid flow approach 
[136' '1371 1138} 11391 1140] . (An alternate strategy involves the construction of gauge-invariant variables: see Refs. 
[141j 142i for reviews.) 

Consider a scalar field in an arbitrary background g^jy. The stress-energy tensor of the scalar field may be written 



(A2) 



Note that we have not yet made any assumptions about the metric or about the scalar field (j). Equation ( A2 ) is 
a completely general expression. We can define a fluid four-velocity for the scalar field by 

(A3) 



It is not immediately obvious why this should be considered a four- velocity. Consider any perfect fluid filling spacetime. 
Each element of the fluid has four-velocity {x) at every point in spacetime which is everywhere timelike, 

w'' (x) {x) = 1 Vx. (A4) 

Such a collection of four-vectors is called a timelike congruence. We can draw the congruence defined by the fluid 



four- velocity as a set of flow lines in spacetime (Fig. 23 1. Each event P in spacetime has one and only one flow line 



passing through it. The fluid four-velocity is then a set of unit-normalized tangent vectors to the flow lines, u^u^ = 1 



For a scalar field, we construct a timelike congruence by Eq. (A3 1, which is by construction unit normalized 



We then define the "time" derivative of any scalar quantity f{x) by the projection of the derivative along the fluid 
four-velocity: 

/ = (A6) 
In particular, the time derivative of the scalar field itself is 



= u^^<j,^^ - ^ff"''</',o'/',/3. (A7) 
Note that in the homogeneous case, we recover the usual time derivative. 



V0 = O^0= ^/(;OV,o</>,o= ^- (A8) 
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\ 



\ 



FIG. 23: A timelike congruence in spacetime. Each event P is intersected by exactly one world line in the congruence. 



The stress-energy tensor (|A2| in terms of (j) takes the form 

'1 



T, 



(A9) 



We can then define a generalized density p and and pressure p by 

1 



P 
P 



1 



(AlO) 



Note that despite the famihar form of these expressions, they are defined without any assumption of homogeneity of 
the scalar field or even the imposition of a particular metric. 



In terms of the generalized density and pressure, the stress-energy (A2) is 
where the tensor /i^jy is defined as: 



(All) 



(A12) 



The tensor h^j^y can be easily seen to be a projection operator onto hypersurfaces orthogonal to the four-velocity 
For any vector field A^^, the product h^^A^ is identically orthogonal to the four- velocity: 



{h^.A'') = A-" (h^^u^) = 0. 



(A13) 



Therefore, as in the case of the time derivative, we can define gradients by projecting the derivative onto surfaces 
orthogonal to the four- velocity 



(A14) 
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FIG. 24: A comoving foliation of spacetime. Spatial hypersurfaces are everywhere orthogonal to the fluid four-velocity it'' 



In the case of a scalar field fluid with four- velocity given by Eq. (A3 1, the gradient of the fleld identically vanishes, 

{Vcjjf = 0. (A15) 

Note that despite its relation to a "spatial" gradient, V/ is a covariant quantity, i.e. a four- vector. 

Our fully covariant definitions of "time" derivatives and "spatial" gradients suggest a natural foliation of the 
spacetime into spacelike hypersurfaces, with time coordinate orthogonal to those hypersurfaces. We can define spatial 
hypersurfaces to be everywhere orthogonal to the fluid flow (Fig. 24 1. This is equivalent to choosing a coordinate 
system for which ~ everywhere. Such a gauge choice is called comoving gauge. In the case of a scalar field, we 
can equivalently define comoving gauge as a coordinate system in which spatial gradients of the scalar field are 
defined to vanish. Therefore the time derivative (A6l is just the derivative with respect to the coordinate time in 
comoving gauge 



d(l> 



(A16) 



Similarly, the generalized density and pressure ( AlO I are just defined to be those quantities as measured in comoving 
gauge. 

The equations of motion for the fluid can be derived from stress-energy conservation, 

T^"',, ^Q^pu^' + [VpY + {p + p) {ii^ + u^'O) , (A17) 
where the quantity O is defined as the divergence of the four-velocity, 

e = (A18) 

We can group the terms multiplied by u'^ separately, resulting in familiar-looking equations for the generalized density 
and pressure 



p + e(p + p) = 0, 

{Wp)>' + {p + p) u'' = 0. 



(A19) 



The first of these equations, similar to the usual continuity equation in the homogeneous case, can be rewritten using 
the definitions of the generalized density and pressure (AlO I in terms of the field as 



4> + + v {(()) = 0. 



(A20) 



This suggests identifying the divergence as a generalization of the Hubble parameter H in the homogeneous case. In 
fact, if we take 5^1/ to be a fiat Friedmann-Robcrtson- Walker (FRW) metric and take comoving gauge, u'^ — (1, 0, 0, 0), 
we have 



3H, 



(A21) 
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and the generalized equation of motion (A20) becomes the familiar equation of motion for a homogeneous scalar, 



+ 3H(j) + V' (0) = 0. (A22) 
Now consider perturbations 5g^„ about a flat FRW metric, 

9t,v = (t) [Vt^u + ^9t,v] , (A23) 

where r is the conformal time and rj is the Minkowski metric rj = diag (1,-1,— 1,-1). A general metric perturbation 
Sgfj.L/ can be separated into components which transform independently under coordinate transformations |jl34j . 

Sg^. = Sgll^'^'- + SglT^^ + SglT""'. (A24) 

The tensor component is just the transverse-traceless gravitational wave perturbation, discussed in Section [V] and 
vector perturbations are not sourced by single-field inflation. We therefore specialize to the case of scalar perturbations, 
for which the metric perturbations can be written generally in terms of four scalar functions of space and time A, B, 
TZ, and Ht- 

Sgao = 2A 
Sgat ^ diB 

Sg,,^2[JZ5,j+d,d,HT]. (A25) 

We are interested in calculating TZ. Recall that in the Newtonian limit of General Relativity, we can write perturbations 
about the Minkowski metric in terms of the Newtonian potential $ as: 

ds^ = (1 + 2$) dt^ - (1 - 2$) S^jdx'dx^. (A26) 

Similarly, we can write Newtonian perturbations about a flat FRW metric as 

ds^ = (t) [(1 + 2$) dr^ - (1 - 2$) Sijdx'dx^ . (A27) 

We therefore expect $ cx 7?. in the Newtonian limit. A careful calculation [ 1381 1141] gives 

a> = -^|^7^, (A28) 



so that in a matter-dominated universe. 



$ = -^7^. (A29) 





In these expressions, TZ is the curvature perturbation measured on comoving hypersurfaces. To see qualitatively why 
comoving gauge corresponds correctly to the Newtonian limit in the current universe, consider the end of inflation. 
Since inflation ends at a particular field value (j) = (pe, comoving gauge corresponds to a foliation for which inflation 
ends at constant time at all points in space: all observers synchronize their clocks to r = at the end of inflation. 
This means that the background, or unperturbed universe is exactly the homogeneous case diagrammed in Fig. [13] 
and the comoving curvature perturbation TZ is the Newtonian potential measured relative to that background. 

To calculate TZ, we start by calculating the four- velocity in terms of the perturbed metric. ^'^ If we specialize to 
comoving gauge, = 0, the norm of the four- velocity can be written 

uf'u^^ a^{l + 2A){u^y = 1, (A30) 
and the timelike component of the four-velocity is, to linear order, 

u'' = -{l-A) 

UQ^a{l + A). (A31) 



This treatment closely follows that of Sasaki and Stewart |139| . except that we use the opposite sign convention for N. 
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The velocity divergence Q is then 



e 



3H 



I- A 



1 fdn l^^dHT\ 



aH \ dr 



dr 



where the unperturbed Hubble parameter is defined as 



Fourier expanding Ht, 



(A32) 



(A33) 



(A34) 



we see that for long- wavelength modes k ^ aH , the last term in Eq. ( A32 1 can be ignored, and the velocity divergence 
is 



I- A 



1 dll 



aH dr 

Remembering the definition of the number of e-folds in the unperturbed case, 



iV 



Hdt. 



(A35) 



(A36) 



we can define a generalized number of e-folds as the integral of the velocity divergence along comoving world lines: 



i j Qds^~^- J e[a{l + A) dr]. 



Using Eq. (A35I for O and evaluating to linear order in the metric perturbation results in 

Af = n- J Hdt, 

and we have a simple expression for the curvature perturbation, 

n = M -N. 



(A37) 



(A38) 



(A39) 



This requires a little physical interpretation: we defined comoving hypersurfaces such that the field has no spatial 
variation. 



(V0)^ = =^ = const. 



(A40) 



Then Af is the number of e-folds measured on comoving hypersurfaces. But we can equivalently foliate the spacetime 
such that spatial hypersurfaces are fiat, and the field exhibits spatial fiuctuations: 



A = n = ^ (j)^ const. 



(A41) 



On fiat hypersurfaces, the field varies, but the curvature does not, so that the metric on these hypersurfaces is exactly 
of the FRW form ( 11 ) with fc = 0. We then see immediately that 



N 



Hdt — const. 



(A42) 



is the number of e-folds measured on fiat hypersurfaces, and has no spatial variation. The curvature perturbation TZ 
is the difference in the number of e-folds between the two sets of hypersurfaces (Fig. 251. This can be expressed to 
linear order in terms of the field variation 54> on fiat hypersurfaces as 

n^M - N = ~^S(I) (A43) 

00 

where TZ is measured on comoving hypersurfaces, and 5N/5<j) and Scf) are measured on fiat hypersurfaces. We can 
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FIG. 25: Flat and comoving hypersurfaces. 



express TV as a function of the field (p: 



N 



For monotonic field evolution, we can express ^ as a function of 0, so that 

SN H 



(A44) 



(A45) 



and the curvature perturbation is given by 



n = j\f-N = - 



SN ^, H 



(A46) 



Note that this is an expression for the metric perturbation TZ on comoving hypersurfaces, calculated in terms of quan- 
tities defined on flat hypersurfaces. For 5(j} produced by quantum fluctuations in inflation, the two-point correlation 
function is 



H 

2^' 



and the two-point correlation function for curvature perturbations is 



277^ TOpiVTre 



(A47) 



(A48) 



which is the needed result. 
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